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FOREWORD 
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A brief  unpublished  report  on  the  subject  of  this  AGARDograph  was  made  by 
Mr  B.Hung(BAe)  at  the  Euromech  75  Coloquium  in  May  of  1976.  The  authors,  who  had 
prepared  that  report,  were  commissioned  by  the  Fluid  Dynamics  Panel  of  AGARD  to  prepare 
the  work  for  this  publication.  The  AGARDograph  Editor  was  Dr  G.G.Pope  (RAE). 
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A COMPARISON  OP  PANEL  METHODS  FOR 
SUBSONIC  PLOW  COMPUTATIONS 
by 

H.A.  Sytsma* 

B.L.  Hewitt** 

P.E.  Rubbert*** 


SUMMARY 

Surface  singularity  or  panel  methods  have,  in 
recent  years,  been  developed  to  a stage  where,  in 
principle  at  least,  they  are  capable  of  providing 
nominally  exact  numerical  solutions  for  incompressible 
potential  flow  around  complicated,  real  aircraft  con- 
figurations. As  such  they  have  proved  to  be  very  useful, 
particularly  to  the  wing  designer. 

There  is  at  present,  a large  variety  of  surface 
singularity  methods  in  use  or  under  development  through- 
out the  world.  In  general,  each  method  has  its  own  model- 
ling, accuracy  limitations,  convergence  characteristics, 
computational  time  per  case  etc..  The  variety  of  methods 
together  with  the  importance  attached  to  them  by  the  wing 
designers  et.al.  means  that  there  is  a real  need  for  a 
data  base  against  which  the  various  programs  (either  exist- 
ing or  under  development)  may  be  checked. 

This  report  contains  such  a data  base  for  a number  of 
relatively  simple  wing  configurations  and  nacelle  configu- 
rations. The  datum  results  have  been  obtained  from  the 
Roberts  (BAe)  Spline-Neumann  Program,  and  a pilot  version  of 
the  Boeing  Advanced  Panel-Type  Influence  Coefficient  Method. 

In  addition,  results  from  the  practical,  engineering 
type  application  of  several  methods  are  compared  with  the 
datum  solutions.  These  comparisons  suggest  that  of  the  methods 
considered,  i.e.  the  NLR  Panel  Method,  Hunt-Semple  "sheets" 
method,  Roberts  Spline-Neumann  program  and  the  pilot  version 
of  the  Boeing  Advanced  Panel-Type  Influence  Coefficient  Method, 
the  latter  is  the  most  efficient,  in  terms  of  the  accuracy/ 
computation  time  ratio,  for  the  simple  test  cases  considered. 
However,  it  must  be  realised  that  these  test  cases  are  not 
representative  of  the  production  cases  which  are  normally 
required  by  aircraft  designers,  and  that  the  calculations 
related  to  the  results  in  this  document  were  carried  out  in 
1976}  many  improvements  to  programs  have  been  incorporated 
since  then. 
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***)  Engineering  Manager,  Boeing  Military  Airplane  Development,  P.O.Box  3999,  Seattle,  Washington  98124,  USA 
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LIST  OF  SYMBOLS 


C local  chord  a free  stream  Incidence  (degrees) 

CLL  1 ocal  lift  coefficient  P sectional  circulation 

Cp  pressure  coefficient  p doublet  strength 

n unit  normal  vector  (outward  positive)  q fractional  spanwise  distance 

0 -r  free  stream  velocity  vector  (J  source  strength 

X,Y,Z  rectangular  co-ordinate  system  0 circumferential  angle  (degrees) 

Vx,V  velocity  components  in  Xf  Y,  Z system  0 perturbation  velocity  potential 

t/c  thiokn  •as/chord  ratio  (P  total  velocity  potential 

C/Dp  chord/exit  diameter  ratio 

1.  INTRODUCTION 

Surface  singular! ty  or  so-called  panel  methods  for  the  computation  of  incompressible  flows  with 
In;  ir  compressibility  rr<  tionn,  have  been  in  use  for  more  than  lh  years.  During  these  years  they  have 
b»  en  levelled  (mainly  because  of  the  advent  of  fast  computers)  into  an  important  and  frequently  used  tool 
for  aerodynamic  i ts. 

Although:  the  lasii-  theorem  (i.c*.  Green’s  theorem)  underlying  all  panel  methods  is  the  same,  there 
exists  at  present  a large  variety  of  panel  methods.  This  is  due  to  several  factors.  First  of  all,  Green’s 
theorem  leave  , in  principle,  a freedom  of  choice  from  a variety  of  combinations  of  singularity  (source, 
loiblet)  list ri but  ions • Secondly,  panel  methods  are  numerical  by  nature,  leaving  for  each  method  devel- 
oper the  possibility  of  choosing  hie  own  numerical  scheme  with  regard  to  e.g.  the  discretisation  of  the 
singularity  distributions,  geometry  representation,  type  of  boundary  conditions  used,  solution  method  etc. 

This  multichoice  situation  naturally  given  rise  to  questions  about  which  type  of  method  is  most  effi- 
cient in  terms  of  the  accuracy/computation  time  ratio  for  a given  type  of  flow  problem.  The  present  report 
present,  the  results  of  a joint  NLR-BAe-Boeing  effort  (with  some  participation  by 
McDonnell  - Douglas.)  to  establish  and  compare  this  "efficiency”  for  a limited  number 
f panel  methods  that  have  been  generally  accepted  in  literature  as  viable  programs  for  aeronautical 
appli  at ions.  Three  of  the  methods  involved  in  the  comparison  are  of  the  "first  - order" 
type  (truncation  error  of  order  (panel  dimension)).  These  methods  are  known  as  the  NLR-Panel  method 
(Ref.  l),  the  Hunt-Semple  "lines  " method  (Ref.  2)  and  the  Hunt-Sample  "sheets"  method  (Ref,  2)  the  last 
two  being  developed  at  BAe  (Warton).  Results  of  a "pseudo-second-order"  method  are  available  from  Boeing, 
using  their  Interim  Higher-Order  method  (developed  by  P.K.  Rubber!  et . al  .),  which  is  a pilot  program  for 
Boeings*  Advanced  Panel-Type  Influence  Coefficient  Method  (Ref,  j). 

The  aim  of  establishing  and  comparing  the  efficiency  of  the  various  numerical  methods  encounters  two 
main  problems.. 

Firstly,  no  exaet.  solutions  currently  exist  for  three-dimensional  lifting  potential  flow  problems. 
However,  in  order  to  obtain  some  kind  of  measure  with  regard  to  accuracy,  the  following  approach  was 
adopted.  For  a selected  set  of,  relatively,  simple  geometries  (a  family  of  wings  and  nacelles)  rather 
" solutions  were  obtained  by  moans  of  the  " third-order"  panel  method  developed  by 

A.  Roberts  at  BAe  (Weybridge)  (fief.  4).  In  this  report  it  is  assumed  that  this  method  provides  results 
with  latum  accuracy  against  which  the  results  of  the  other  methods  can  be  checked.  Results  with  datum 
accuracy  were  also  provided  by  Boeing. 

.V  ondly,  there  remains  the  problem  of  comparing  computing  times.  Computing  costs  related  to  the 
methods  are  very  difficult  to  compare  in  a truly  fair  manner.  Even  if  it  were  possible  to  run  all  pro- 
grams individually  on  the  same  computer  with  access  to  the  full  storage  capacity,  the  computation  times 
•ould  only  give  a rather  narrow  comparison.  For  example,  a versatile  program  would  probably  use  more  com- 
putation time  than  a specialised  one.  In  the  comparative  study  being  described  the  situation  was  even  more 
complicated,  since  the  calculations  were  carried  out  on  different  computers.  Nevertheless,  a comparison  of 
’imputation  times  is  presented.  This  was  constructed  tyr  scaling  all  computation  times  to  a reference  com- 
puter, which  in  this  case  wan  chosen  as  the  CDC  6600.  However,  it  will  be  clear  that  the  scaling  factors 
should  be  interpreted  with  care,  and  therefore  the  actual  computation  times  for  each  computer  are  also 

•)  Withii  the  mope  of  this  report,  converged  should  be  interpreted  in  the  sense  that,  an  increase  of  the 
number  of  panels  used  did  not  affect,  the  results  significantly. 
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presented. 

In  the  next  chapter  the  mathematical  background  of  the  calculation  of  inviscid  incompressible  flow 
around  arbitrary  lifting  bodies  will  be  briefly  outlined.  Some  details  about  the  particular  numerical 
schemes  used  in  the  methods  being  currently  compered  are  presented  in  chapter  3.  The  geometrical  details 
of  the  selected  set  of  test  configurations,  together  with  flow  conditions  and  the  chosen  panel  distribu- 
tions are  given  in  chapter  4.  A discussion  of  the  results  in  terms  of  chordwise  pressure  distributions, 
velocity  component^  quantities  such  as  sectional  load  and  circulation  is  presented  ir.  chapter  5*  The  sen- 
sitivity of  the  results  to  the  number  of  panels  used  is  subject  of  chapter' 6.  In  chapter  7 an  attempt  is 
made  to  relate  accuracy,  number  of  panels  used,  and  computation  time  for  the  various  methods  compared. 
Finally  some  concluding  remarks  are  given  in  chapter  8. 

Some  provisional  results  of  this  study  were  presented  by  B.  Hunt  (BAe)  at  the  Euromech  75  colloquium 
at  Braunschweig  in  May  '76.  It  was  felt  by  the  participants  to  be  worthwhile  to  make  the  results  available 
to  a wider  audience,  and  consequently  publication  through  the  AGARD-FDP  was  sought  and  granted. 


2.  THBOHETICAL  BACKGROUND  OF  PANEL  METHODS 

In  this  chapter  the  mathematical  basis  for  the  treatment  of  3-D  irrotational  incompressible  flow 
around  arbitrary  configurations  will  briefly  be  outlined.  Only  this  type  of  flow  is  considered  in  this 
report.  More  details  of  the  theory  can  be  found  in  standard  references  (Refs.  5*6,7)  • 

Flows  of  this  type  are  characterized  by  a perturbation  velooity  potential  9.  This  potential  satisfies 
Laplace's  equation: 

9 + 9 +9  =0  (l) 

in  a region  R surrounding  the  body,  and  is  subject  to  cei*ain  boundary  conditions.  The  body  is  bounded  by 
the  surface  S.  In  order  to  obtain  a unique  solution  in  lifting  cases,  a potential  discontinuity  surface  W 
(the  so-called  wake),  leaving  the  sharp  trailing  edge  and  extending  to  infinity,  has  to  be  introduced 
(Fig.  l). 

According-to  Green's  third  identity,  a solution  of  eqn.  (l)  may  be  expressed  at  any  point  P as  the 
potential  induced  by  a combination  of  so-called  source  singularities  of  strength  O and  doublet  singulari- 
ties  of  strength  n,  distributed  on  the  surfaces  S and  W: 

o(p)  - as  + Jju(Q)  (ji-;)  as  (2) 


where  r is  the  distance  from  the  field  point  P to  the  surface  paint  Q,  and  is  the  derivative  in  the 
direction  of  the  outward  surface  normal.  A required  solution  may  be  found  by  imposing  suitable  boundary 
conditions  on  S and  W.  Within  the  scope  of  this  report  only  the  Neumann  type  of  boundary  condition  is  con- 
sidered, which  can  be  imposed  directly: 

-Voo  ■ "Q  (3) 

or  (indirectly)  by  requiring  the  Dirichlet  condition  that 

9 » 0 (4) 

on  the  inner  surface  of  the  body.  Upon  making  an  "a  priori"  choice  for  the  doublet  strength  distribution 
when  using  eqn.  (3),  or  an  "a  priori"  choice  for  the  source  strength  distribution  when  using  eqn.  (4)1  a 
Fredholm  integral  equation  of  the  2nd  kind  for  the  remaining  unknown  singularity  strength  is  obtained*). 

It  should  be  noted  that  for  flows  with  circulation,  doublet  singularities  must  be  used  somewhere  in  the 


Common  to  all  panel  methods  is  the  subdivision  of  the  surfaces  S and  W into  so-called  panels  that 
approximate  the  geometric  surfaces  to  a certain  order.  Further  the  singularity  strengths  may  be  chosen  to 
vary  in  a convenient  prescribed  way  over  each  panel,  e.g.  in  a "first-order"  method,  constant  strength 
source/doublet  strengths,  or  possibly  bi-linear  doublet  strength  variations  are  employed  on  flat  panels. 

Selecting  a number  of  control  points  per  panel  (this  number  depending  on  the  order  of  the  singularity 
distribution)  i.e.  points  where  the  boundary  condition  is  applied,  leads  (after  integration)  to  a linear 
system  of  algebraic  equations  in  the  unknown  singularity  strengths.  This  system  may  then  be  munerically 
solved  in  several  ways  e.g.  directly  through  matrix  inversion  or  through  some  iterative  process. 


*)  Other  choices  may  lead  to  an  integral  equation  of  the  1st  kind 
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The  above  brief  description  contains  the  basic  elements  underlying  all  the  panel  methods  involved  in 
the  current  comparative  study*  In  the  next  chapter  these  methods  will  be  described  in  more  detail,  paying 
regard  to  the  schemes  adopted  for  geometry  approximation,  the  location,  type  and  variation  of  singularities 
used,  the  type  of  boundary  conditions  employed  and  the  formation  and  solution  of  the  linear  systems  of  al- 
gebraic equations* 


3.  SHORT  OUTLINE  OP  THE  PANEL  METH0D6  INVOLVED  IN  THE  COMPARISON 
3.1  Introduction 

In  this  chapter  the  methods  involved  in  the  current  comparison  will  be  discused  in  some  detail.  More 
details  about  a particular  method  or  its  scope  of  applications  rnay  be  found  in  the  relevant  references.  It 
should  be  noted  that  the  description  relates  to  the  situation  at  the  time  that  the  calculations  were  rai- 
ned out  (1976).  Since  then  some  of  the  methods  have  been  changed,  or  have  evolved  in  some  respect. 


3 %2  The  NLR  Panel  Method 

The  NLR  Panel  Method  follows,  at  least  with  regard  to  the  solution  of  incompressible  lifting  flow 
problems,  the  lines  of  the  method  developed  by  P.E.  Rubbert  et.  al.  (Ref.  3).  This  method  is  a "first  o — 
derH  method,  and  this  implies  that  generally  the  truncatation  error  is  of  0(h),  where  h is  a characteristic 
panel  dimension.  It  should  be  noted,  however,  that  for  reasons  explained  below  the  error  becomes  of  0(l)  in 
the  limiting  case  of  vanishing  thickness. 

In  this  method  the  body  and/or  wing  surfaces  are  approximated  by  flat  panels  in  the  same  way  as  origi- 
nally defined  by  Hess  and  Smith  (Ref.  )) . 

The  singularities  used  in  the  NLR  Panel  Method  are:  (i)  constant  strength  sources  distributed  on  each 
body  and/ox  wing  surface  panel,  and  (n)  constant  strength  doublet  panels  on  all  wake  surfaces  of  lifting 
on  figuration  8.  For  numerical  reasons  the  wing  wake  is  extended  inside  the  wing,  along  the  camber  surface. 
Prom  the  well-known  induced-velocity  equivalence  between  a constant  strength  doublet  panel  and  a line 
vortex  along  the  edge  of  the  panel,  it  follows  that  the  constant  doublieity  panels  on  wake  and  camber  sur- 
face may  be  re-interpreted  as  a vortex  network  or  lattice.  For  the  doublet  panels  on  the  camber  surface, 
the  chap*-  of  the  chordwise  variation  of  doublet  strength  (which  is  equivalent  to  the  valuation  of  the  in- 
ternal vortex  strengths)  is  prescribed  and  kept  constant  span wise.  However,  an  unknown  scaling  factor  is 
associated  with  each  streamwise  strip  and  is  effectively  determined  by  locally  applying  the  Kutta  condition. 
In  the  NLR-method  a standard  parabolic  shape  is  used  for  the  chordwise  variation  of  discrete  vortex 
strengths  (see  Pig.  2).  As  a precaution  the  suitability  of  this  shape  may  be  checked  for  a typical  4-D  sec- 
tion of  the  wing  under  consideration. 

The  surface  boundary  condition  is  of  the  Neumann  type,  i.e.  (eqn,  (i)) 


dn 


where  n is  the  outward  normal  to  the  body,  and  lly,  the  free  stream  velocity  vector.  This  b.c.  is  applied  at 
the  geometric  centroid  of  each  flat  surface  panelf  i.e.  at  the  so-called  "collocation  points".  In  lifting 
■uses  the  wake  doublet  strength  is  determined  by  forcing  the  flow  to  leave  the  wing  tangential  to  a plane 
that  bisects  the  wing  trailing  edge  at  a small  ( - 10  local  chord)  distance  downstx'eam  (i.e.  the  Kutta- 

condition,  see  Pig. 3).  The  general  source  panel  and  vortex  lattice  arrangement  is  depicted  in  figure  4. 

The  matrix  equation  stating  that  the  unknown  source  and  doublet  strengths  must  satisfy  the  Neumann 
boundary  condition  ( eqn • (3 ) ) at  all  collocation  points,  together  with  the  Kutta— condi tion,  may  be  written 
in  matrix  form  as  follows: 


[a]  [x]  . [b] 

J 7 x alj  XJ  - bi  (i-1»N) 


(M 

(6) 


where  the  are  so-called  normal  velocity  influence  coefficients,  expressing  the  normal  velocity  at  col- 
location point  i due  to  a unit  strength  singularity  on  panel  j.  Here  stands  for  the  boundary  condition 
pertaining  to  collocation  point  i,  and  x.  for  the  unknown  singularity  strength  on  panel  j*  The  normal  velo- 
city influence  coefficients  are  calculated  and  stored  on  disc  column-wise,  since  then  only  the  data  per- 
taining to  one  panel  has  to  reside  in  central  memory.  The  velocity  components  V^,  and  V%  are  likewise 
calculated  and  stored  columnwise.  Use  is  made  of  fax—  and  neat— field  approximations  given  by  Hess  and  Smith 


(Rpf.9). 
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The  linear  system  of  algebraic  equations  is  solved  iteratively.  For  this  purpose  the  matrix  £a]  and  the 
vectors  [x]  and  [b]  are  organized  in  a special  way,  as  depicted  in  figure  5.  The  coefficients  relating  to 
a streamwise  source  panel  strip  are  grouped  together  in  blocks  (within  the  sub-matrix  [Sc]  ) on  the  main 
diagonal.  This  structure  allows  the  solution  of  the  system  by  means  of  an  adapted  (note  the  column  wise 
storage  of  the  influence  coefficients)  block  Gauss-Seidel  approach.  Convergence  is  tested  on  Afj  and  Af 
for  two  consecutive  iteration  steps  (A(j  » (J  n+*  - £7n,  n being  the  iteration  number).  The  process  is  stop- 
ped when  both  A O and  AF  are  < .0001,  For  a simple  lifting  wing  configuration  a number  of  10  -15  itera- 
tions is  typical. 

More  details  about  the  method  with  regard  to  the  formation  and  the  solution  of  the  system  of  equations 
can  be  found  in  reference  10,  Reference  11  contains  a program  listing  on  micro-fiche.  Examples  of  applica- 
tions can  be  found  in  references  12,  13, 

3*3  Hunt  - Semple  Panel  Methods 

The  Hunt  - Semple  "Mark  1A"  panel  program  contains  within  a single  program  the  option  of  using  a 
"lines"  model  or  a "sheets"  model. 


3.3.1  The  "lines"  model 

The  Hunt  - Semple  "lines"  method,  like  the  NLR-Panel  method,  is  a "first-order"  method.  Both  the  NLR- 
method  and  the  Hunt  - Semple  "lines"  method  use  the  same  type  of  singularity  distributions  and  boundary 
conditions  (see  Fig.  4),  but  differ  in  two  important  respects:  (i)  in  the  method  of  applying  the  Kutta- 
condition  and  (ii)  in  the  Hunt  - Semple  use  of  an  "optimiser"  to  determine  an  "optimal"  chordwise  shape  for 
the  strengths  of  the  internal  constant  doublicity  panelsf  n.b.  a new  shape  is  calculated  for  each  wing. 

In  the  Hunt  - Semple  "lines"  panel  method  the  following  approach  is  adopted  for  determining  a charac- 
teristic "optimal"  variation  or  shape  of  the  chordwise  doublet  strength.  It  is  argued  that  for  a three- 
dimensional  lifting  component  this  chordwise  shape  may  be  conveniently  chosen  from  a quasi-two  dimensional 
argument.  For  each  lifting  strip  the  set  of  doublet  strengths  is  chosen  such  that  for  an  onset  flow  of 
90°,  with  the  internal  vortices  extended  indefinitely  spanwise,  and  all  other  influences,  including  the 
source  panels  neglected,  the  doublet  distribution  on  the  camber  line  satisfies  (in  a weighted  least  squares 
sense),  the  surface  boundary  conditions  at  the  surface  collocation  points.  An  example  of  such  an  "optimised" 
distribution  is  given  in  figure  6a.  The  corresponding  vortex  distribution  (Fig.  6b)  can  differ,  signifi- 
cantly, from  the  type  used  at  NLR  (see  Fig.  2),  and  from  the  distribution  usually  used  at  MBB  (Ref.  14) 
where  the  strengths  are  chosen  proportional  to  local  thickness. 

In  the  currently  described  Hunt  - Semple  methods,  in  order  to  make  use  of  a 2-D  "optimised"  chordwise 
doublicity  shape  1:  3”D  calculations,  the  shape  values  of  doublicity  along  each  chordwise  strip  are  scaled 
by  the  associated  trailing  edge  value  of  doublicity^  which  is  regarded  as  the  unknown  for  each  chordwise 
strip.  Figure  6c  shows  the  source  strengths  associated  with  these  "optimised"  doublet  distributions.  It 
can  be  seen  that  the  "optimised"  distribution  produces  much  smaller  source  strengths  than  e.g.  does  the 
thickness— based  distribution.  The  objective  of  this  procedure  is  to  reduce  the  discretisation  errors  due 
to  the  fairly  crude  constant  strength  source  panel  modelling.  It  is  to  be  noted  that,  thinking  similar  to 
that  underlying  the  Hunt  - Semple  "optimiser"  is  expressed  by  Rubbert  et.  al.  in  Reference  15  • 

With  respect  to  satisfying  the  Kutta-condition  the  following  remarks  can  be  made.  In  the  Hunt  - Semple 
"lines"  method  a set  of  Kutta-points  is  first  constructed  downstream  of  the  wing  trailing  edge  on  the  ex- 
tension of  its  camber  surface.  Whereas  in  the  NLR-method  a tangential  flow  condition  is  fulfilled  at  such 
points,  the  technique  adopted  here  could  be  described  as  a "calculated"  Kutta-condition. 

The  approach  is  based  on  an  exact  analysis  of  the  2-D  flow  around  a particular  family  of  analytical 
airfoils,  and  is  characterized  by  the  chord  length  C and  the  trailing-edge  angle  A.  The  following  function 
has  been  derived  describing  the  velocity  component  normal  to  the  trailing-edge  bisector  as  a function  of 

A,  the  circulation  F and  the  distance  6 downstream  of  the  trailing  edge. 

T + A 


■ "tc  <§>2*-  A 


(7) 


This  formula  is  applied  directly  in  the  3-D  program,  where  the  Kutta-points  are  placed  at  about  ^ 


of  the 

local  chord  downstream  of  the  trailing  edge  along  the  mean  spanwise  position  of  the  surface  collocation 
points  in  the  relevant  atrip.  The  reasoning  behind  this  approach  is  to  allow  the  Kutta-points  to  be  suffi- 
ciently far  removed  from  the  trailing  edge  that  flow  field  induced  at  these  points  by  the  actual  constant 
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at ren^th  source  panels  closely  matches  that  whic)  would  be  induced  by  a corresponding  "higher-order"  sour- 
ce distribution.  The  circulation  used  in  eqn  (?)  is  found  directly  as  part  of  the  solution  of  the  system 
of  linear  rquations  represented  by  eqn  (6),  The  solution  of  this  system  depends  on  the  nonnal  velocity  in- 
duced by  all  the  singularities  at  collocation  points.  Thus  because  of  implied  small  source  strengtns  the 
normal  velocity  induced  by  the  internal  vortex  system  is  of  primary  importance.  Now  it  car  be  rlxnti  flftef.8)  that 
the  normal  velocity  induced  by  a discrete  vortex  system  is  much  less  sensitive  to  the  parameter  h/A  than 
e.g.  the  tangential  velocity  (Pig.7).  Fi  om  this  it  also  follows  that  the  circulation  I'  will  be  less  sensi- 
tive to  a decrease  in  h/A  than  the  sectional  integrated  load.  More  details  about  this  approach  can  be 
found  in  reference  2. 

The  matrix  formation  is  essentially  the  same  as  in  the  NLR-method,  except  that  the  matrix  is  written 
to  disc  row-wi se,  and  the  chordwise  chain-  of  the  doublicity  variation  may  vary  spanwise  (n.b.  usually  this 
: c-  mill  not  uiwi  u.-l  is  calculated  by  the  so-called  "optimi ser",  a • irlier). 

fht  linear  aquations  a>-»*  solved  by  a point  Gauss— Seidel  nethodi (or optionally,  by  over- or  underrelaxation) 
except  that  the  influences  of  the  vortex  chordwise  vortex  "ladders**  on  the  Kutta-points  are  inverted  as  a 
single  diagonal  block;  as  in  the  MBB-method  (see  Ref.  14).  The  solution  is  considered  to  have  converged 
when  the  residues  (b^  - / a(  x^.)  computed  during  an  iterative  cycle  are  all  less  than  .0003. 


3*3»2  The  "sheets"  model 

The  Hunt  - Semple  "sheets"  method  differs  from  the  "lines"  method  only  in  the  sense  that,  on  the 
(flat)  camber  surface  panels,  the  doublicity  on  each  panel  is  built  up  from  a linear  combination  of  four 
bilinear  Lagrangian  interpolation  modes,  one  mode  being  associated  with  each  of  the  four  panel  comers.  A 
main  requirement  for  the  currently  described  "sheets"  model  is  to  provide  sensible  estimates  for  both  the 
"tangential"  first  derivatives  of  doublicity,  and  the  doublicity  values,  at  panel  centroids  (see  Fig.  8). 

For  the  purpose  approximating  the  tangential  derivatives,  doublicity  values  at  adjacent  centroids  are  lo- 
cally fitted,  three  at  a time,  by  quadratics.  Independent  fits  are  made  in  both  the  chordwise  and  spanwise 
directions.  These  local  fits  provide  analytical  estimates  for  first  derivatives  of  doublicity  at  panel 
centroids,  which  involve  doublicity  -allies  at  adja  ent  centroids.  An  described  previously,  the  only  un- 
known associated  with  the  doublicity  valuation  al  mg  each  chordwise  strip  of  panels  is  the  value  at  the 
trailing  edge.  Thus,  using  the  first  deriv. estimates  there  results  a doublicity  representation  over 
each  panel  which  depends,  in  genera'.,  on  three  unknown  trailing  edge  values  of  doublicity.  This  doublicity 
is  equivalent  to  distributed  vorticity  - n x grad  u on  each  panel,  together  with  a discrete  line  vortex 
of  linearly  varying  strength  u along  each  panel  edge.  It  can  be  shown  that  for  the  "sheet"  approach,  par- 
ticularly in  the  case  of  lifting  wings  with  small  thickness,  the  tangential  velocity  at  surface  collocation 
points  will  be  better  represented  than  for  the  "lines"  model,  where  a discrete  internal  vortex  lattice  is 
employe!  . The  effect  is  demonstrated  for  a simplified  model  in  figure  7*  Further  it  should  be  noted  that 
due  to  the  discretisation  of  the  lifting  vortex  system  in  a "lines"  model,  the  jump  in  the  spanwise  velo- 
ity  component  at  the  trailing  edge  can  never  be  predicted  correctly,  whereas  a "sheets"  model  allows  for 
this  discontinuity.  With  regard  to  the  errors  for  a "lines"  model  associated  with  the  chordwise  and  the 
spanwise  velocity  components  at  the  wing  surface,  it  should  be  noted  that  the  parameter  A from  figure  7 
should  be  interpreted  respectively  as  the  chordwise  panel  length  (error  in  V^),  or  as  the  spanwise  strip 
width  (error  in  V ).  It  may  be  noted  that  these  errors  in  the  velocity  mav  lead  to  an  error  in  the  pres- 

y 

sure  distribution  and,  indirectly,  to  an  error  in  integrated  quantities  such  as  the  sectional  load.  This 
does  not  apply,  however,  to  the  circulation,  as  explained  earlier  in  3*3*1» 

The  chordwise  share  of  the  doublicity  on  the  camber  surface  panels  is  found  through  essentially  the 
same  "optimising"  process  as  for  the  "lines"  model.  It  should  be  noted,  however,  that  because  of  the  2-D 
nature  of  the  optimisation  process,  the  truncation  error  for  vanishing  thickness  in  the  3—D  case  is  still 
formally  of  0(l),  as  for  other  "first-order"  panel  methods. 

The  Kutta— condition  is  the  same  "calculated"  one  as  employed  in  the  "lines"  method,  and  was  outlined 
in  chapter  3*3.1 • 

The  matrix  formation  for  the"sheets"  model  is  somewhat  more  complex  than  that  for  the  "lines"  model. 

As  previously  explained,  the  doublicity  on  each  panel  is  generally  dependent  on  three  unknown  trailing  edge 
values  of  doublicity.  This  dependence  carries  through  directly  into  the  panel  influence  expressions,  and 
it  is  this  fact  which  slightly  complicates  the  "sheets"  model  matrix  formation  procedure  relative  to  that 
for  the  "lines"  model;  further  details  are  given  in  reference  7. 
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The  solution  of  the  linear  system  of  equations  le  the  same  as  for  the  "lines"  model, 

3,4  Boeing’s  Interim  "Higher-Order"  Panel  Method 

The  type  of  geometry  approximation  in  Boeing’s  method  is  "first-order",  i,e.  flat  panels  were  employed 
or.  the  surface  and  the  wake  of  the  test  configurations.  The  wake  is  not  extended  inside  the  wing.  Generally, 
a configuration  is  divided  into  so-called  "networks".  A "network"  is  viewed  as  a portion  of  the  boundary 
surface  which  is  subdivided  into  panels,  and  is  logically  independent  in  the  sense  that  it  contributes  as 
many  equations  to  the  overall  problem  as  it  contributes  unknowns. 

In  contrast  to  the  methods  discussed  in  the  preceding  sections,  for  this  method  the  singularity  dis- 
tribution on  the  surface  (and  wake)  panels  is  of  doublet  type  only.  The  doublet  strength  distribution  on 
each  surface  panel  is  chosen  to  vary  quadrati rally  in  two  orthogonal  directions^  a local  2-D  Taylor  expan- 
sion involving  6 parameters.  Discirte  values  of  the  doublet  strength  are  assigned  to  certain  standard  points 
on  each  "network".  The  location  of  these  points  (l  ,e,  doublet  parameter  points)  is  shown  in  figure  9*  The 
doublet  distribution  on  a surface  panel  is  then  found  by  fitting  the  6 parameter  quadratic  form,  in  a 
i gh t ed  least-squares  sense,  using  the  doublioity  value  at  its  centre  point  and  at  those  of  the  adjacent 
panels  ( ) points  in  all).  The  weight  is  chosen  to  be  relatively  very  large  only  for  the  doublet  value  point 
on  the  panel  considered.  For  a panel  adjacent  to  a "network"  edge,  doublet  value  points  on  the  edge  are  al- 
so included  in  the  fit  (Fig.  9). 

The  control  points  on  a network,  i.e.  points  where  boundary  conditions  are  applied,  are  also  indicated 
in  figure  9*  These  points  include  panel  centre  points  as  well  as  edge  points  (very  slightly  displaced  from 
"r.  e t wo  rk  " edges  ) • 

At  panel  centre  points  a Dirichlet  condition  is  employed  which  requires  the  vanishing  of  the  total 
potential  on  the  interior  surface  of  the  wing.  The  edge  control  points  serve  to  facilitate  matching  between 
"networks".  Thus  at  network  edges,  continuity  of  the  doublet  strength  and  its  gradient  is  ensured  to  a cei— 
tain  order.  At  the  trailing  edge  (where  the  wing  and  wake  "networks"  abut)  the  Kutta-condition  is  satis- 
fied implicitly.  For  such  a model  the  local  doublet  strength  ma y be  identified  with  the  total  potential <D  $ 
which  is  the  basic  unknown  in  this  particular  version  of  the  Boeing’s  method. 

This  approach  leads  to  a linear  system  of  algebraic  equations  in  the  unknown  doublet  strength  para- 
meters. This  system  is  solved  by  employing  a Crout  decomposition  algorithm  with  pivoting  in  diagonal  blocks 
only.  Also  use  is  made  of  a very  fast  Compass-coded  vector  product  subroutine.  After  solving  this  system 
the  velocity  components  at  any  point  are  directly  found  from  the  analytic  gradient  of  the  total  potential 

<p . 

2 

Note  that  the  truncation  error  of  the  approach  described  above  is  generally  0(h),  but  becomes  0(h  ) 
in  case  of  vanishing  surface  curvature  (h  being  a characteristic  panel  dimension). 

More  details  about  this  method  can  be  found  in  reference  3* 
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3.5  Roberts’  Spline  - Neumann  Method 

In  Roberts’  method  (Ref.  4)  the  aircraft’s  wetted  surfaces  are  subdivided  into  a set  of  convenient 
"carpets",  which  are  analogous  to  the  "networks"  in  the  Boeing  method  (see  chapter  3»4).  Each  carpet  is 
subdivided  into  curvilinear  panels  whose  comers  form  a grid  of  points.  These  carpets  are  first  mapped  in- 
to a set  of  rectangles  in  convenient  parametric  u,  v planes.  The  Cartesian  components  of  the  vector  posi- 
tion of  each  grid  point  are  specified  as  part  of  the  input.  Each  of  these  components  is  regarded  as  a re- 
gular bi-cubic  spline  function  of  the  parameters  u and  v in  the  parametric  plane.  Using  this  representation, 
the  surface  shape  associated  with  any  single  panel  is  defined  by  sixteen  basic  bi-cubic  spline  modes  (one 
mode  is  shown  in  Fig.  10).  Thus  continuity  upto  the  second  derivative  between  panel  edges  is  preserved. 

The  same  approximation  is  used  to  represent  the  rigid  wake  surface  and  its  extension  inside  the  wing. 

The  type  and  location  of  the  singularities  used  in  Roberts*  method  correspond  to  those  of  the  "first- 
order"  NLR  and  Hunt  - Semple  methods.  Sources  are  employed  on  the  wetted  surfaces  and  doublioity  on  the 
wake  and  its  extension  along  the  camber  surface  inside  the  wing.  As  in  the  Hunt  - Semple  and  NLR  methods 
the  ohordwise  shape  of  the  internal  doublet  distribution  is  prescribed.  Roberts’  chordwlse  shape  is  given 
by  the  integral  of  a mixture  of  Bimbaum  vorticity  modes.  The  unknown  source  and  doublet  strength  distri- 
butions over  the  set  of  "carpet"  rectangles  in  the  parametric  plane  are  cubic.  Again  a bi  -cubic  spline 
approach  is  used  and  any  single  panel  is  covered  by  sixteen  basic  spline  modes  (see  Fig.  10).  By  applying 
the  boundary  condition  of  zero  normal  velocity  at  the  comers  of  each  panel,  a one-to-one  correspondence 
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with  the  unknown  singularity  strength  in  the  middle  of  each  mode  (point  A in  Fig.  10 ) in  obtained. 

In  contrast  with  the  "fi rst— order”  methods  described  earlier  the  Kutta— condi ti on  (i.e.  the  require- 
ment that  velocities  remain  finite  at  the  wing  trailing  edge)  is  satisfied  implicity.  Over  a certain  num- 
ber of  panels  on  both  the  upper  and  lower  surfaces  of  the  wing  adjacent  to  the  trailing  edge,  special  sin- 
gular types  of  source  strength  modes  are  employed.  These  modes  are  based  on  the  analysis  of  a wedge  type 
of  flow,  thus  implicitly  ensuring  the  correct  flow  behaviour  at  the  trailing  edge.  Wing  tips  are  treated  in 
a similar  way. 

The  influence  coefficient  of  each  singularity  mode  for  the  normal  velocity  at  each  collocation  point 
ia  found  numerically  through  specially  developed  Gauss  quadrature  techniques. 

The  resulting  linear  system  of  algebraic  equations  is  solved  directly  by  Grout  matriA  inversion. 

Note  that  the  truncation  error  of  this  method  is  generally  of  0(h  ) (h  being  a characteristic  panel  dimen- 
sion) but,  because  of  the  fixed  choice  of  the  internal  doublet  distribution,  this  becomes  formally  of  0(l) 
in  the  liriting  case  of  3“P  lifting  wings  with  vanishing  thickness. 

More  details  can  be  found  in  reference  4. 

4.  DEFINITION  OF  THE  TEST  CASES 

4.1  Definition  of  the  Geometry  of  the  Test  Configurations 

Three  different  types  of  relatively  simple  test  configurations  were  chosen  for  the  current  comparison. 
The  first,  configuration  is  a swept  tapered  wing  with  RAE  wing  "A"  planform,  without  camber  or  twist  and 
without  dihedral.  The  second  one  is  a more  complex  configuration  of  the  st raked— wing  typo.  The  third  confi- 
guration considered  is  an  annular  duct  or  flow-through  nacelle.  This  latter  configuration  was  introduced 
because  it  has  been  observed  that  the  computation  of  such  "partially  internal"  flow  poses  accuracy  pro- 
blems to  several  panel  methods.  The  plan  forms  of  these  wings  and  the  crossection  of  the  nacelle  are  defined 
in  figure  11. 

For  all  configurations  the  same  airfoil  section  was  used  for  simplicity  of  definition,  this  was  taken 
from  the  NACA-Four-Digi t series.  This  appeared  to  be  of  importance  for  the  methods  of  Boeing  (Rubbert)  and 
Roberts  since  these  methods  cannot  easily  handle  open  trailing  edeges  without  including  :\  thick  wake  repre- 
sentation. Consequently  the  trailing  edges  had  to  be  handled  in  a special  way.  Rubbert  solved  this  problem 
by  altering  the  slope  of  the  trailing  edge  panels,  whereas  Roberts  extrapolated  the  trailing  edge  slight- 
ly beyond  the  original  trailing  edge  to  enable  closure. 

In  order  to  study  the  capability  of  the  various  methods  to  treat  wings  of  different  thicknesses,  a 
"family"  of  wings  was  defined,  i.e.  calculations  were  carried  out  for  throe  thickness/chord  ratios  for  the 
RAE  WING  planform  (t/c  - .15 » .05  and  .02)  and  for  two  T/C  ration  for  the  st raked— wing  plan forms 
(T/C  - .05  and  .02). 

Although  the  extremely  thin  wing  configurations  (T/C  . ,0?)  are  not  "realistic"  in  the  sense  that 
wings  of  this  thickness  arc  not  practical  from  a constructional  point  of  view,  they  were  introduced  for 
the  following  reason.  In  the  computation  of  flows  with  allowance  for  linear  compressibility,  the  correc- 
tion, according  to  Goethert's  rule,  requires  the  computation  of  the  incompressible  flow  around  an  effecti- 
vely thinner  configuration.  This  means  o.g.  that  for  a 5 % thick  configuration  at  M .8,  the  incompressi- 
ble flow  in  calculated  around  a configuration  only  \ 7*  thick.  Also  two  nacelles  were  studied,  with  two 
chord/exit  diameter  ratios,  i.e.  1,0  .and  1.333. 

4.2  Flow  Conditions 

With  one  exception  the  calculations  were  carried  out.  for  (t  51  and  with  a rigid  wake.  E’or  the  wing 
casen  the  latter  war.  located  in  the  horizontal  plan*'  of  symmetry  downstream  the  wing.  The  wake  of  the  na- 
celles is  located  on  a cylindrical  surface  with  its  diameter  equal  to  the  exit  diameter  of  the  nacelle,  Oi- 
ly one  wing  care  (i.e.  RAE'  WING  , T/o  .15)  in  included  for  <1  0 since  it.  can  be  argued  that  none 

of  the  methods  involved  should  have  problems  al  uliting  non-lifting  potential  flow. 

4.3  Panel  Distributions 

Apart,  from  defining  the  test  con  figurations  geometrically,  it.  war  also  deemed  necessary  to  require 
the  use  of  identical  panel  distributions  for  it  least  the  "first-order"  methods.  Both  the  spanwise  and  the 
chordwisr  panelling  was  therefore  prescribed.  The  two  chosen  distributions  are  tabulated  in  tabic  1 (chord- 
wise)  and  table  2 (spanwise).  In  thin  way,  differences  between  the  results  will  be  revealed  more  clearly 


1 

j 


-9- 


and  will  not  be  obscured  by  the  use  of  different  panel  distributions,  Thus,  some  useful  information  will 
also  be  obtained  about  the  sensitivity  of  the  results  to  the  number  of  panels  used  in  the  computations. 

Panelling  for  the  "higher-order*  methods  was  loft  open  to  choice  in  order  to  create  the  best  conditions 
for  obtaining  supposedly  datum  accuracy.  Furthermore,  it  was  also  requested  that  the  "higher-order"  methods 
provide  so-called  "engineering  solutions"  (i,e,  results  with  acceptable  accuracy  m U;rmu  of  chordwiue  pres- 
sure distributions)  again  with  freedom  of  panelling, 

5,  EVALUATION  OF  THE  HBSULTS 

5.1  Datum  Results 

5.1.1  Introductory  Remarks 

In  this  chapter  datum  solutions  from  the  methods  of  Roberts  and  Rubbert  will  be  discussed.  The  results 
are  compared  graphically  in  terms  of  chordwiae  pressure  distributions,  velocity  distributions,  and  circula- 
tion. Unfortunately  sectional  load  distributions  could  not  be  compared  since  Roberts*  program  did  not  pro- 
vide this  data.  Corresponding  tabulated  data  can  be  found  in  table  3* 

The  panel  distributions  used  to  generate  the  datum  results  were  the  following.  Rubbert  employed  40 
panels  chordwiae  and  12  apanwisc  strips  for  the  RAE  WING  cases.  Roberts  used  39  panels  chordwise  and  13 
apanwise  strips  for  these  cases. 

For  the  STRAKKD  WINO  cases  Rubbert  employed  38  panels  chordwise  and  12  spanwise  strips,  whereas 
Roberts  used  39  panels  chordwise  and  18  spanwise  strips  ( thickness/chord  ratio  .05)  and  24  strips  spanwise 
for  T/C  - .02.  On  the  nacelles,  the  number  of  chordwise  panels  used  were  40  and  55,  by  Rubbert  and  Roberts, 
respectively.  Both  used  10  circumferential  strips. 

It  should  be  noted  that  the  set  of  data  is  not  complete,  in  the  sense  that  datum  results  from  both 
methods  are  not  available  for  all  test  cases.  With  regard  to  Rubbert *s  results  it  should  be  noted  that  an 
anomaly  in  the  immediate  vicinity  of  the  trai ling-edge  is  present;  this  is  due  to  the  alteration  of  the 
trailing  edge  elopes  to  enable  closure  (see  4.1 )•  The  consequences  of  this  procedure  were  checked  two-di- 
mensional ly  and  appeared  to  have  no  other  effects. 

5.1.2  Discussion 
RAE  WINO  Cases 

Results  for  the  RAE  WINO  cases  are  presented  graphically  in  figures  12  to  33.  Chordwise  pressure  dis- 
tributions show  generally  a high  level  of  agreement  for  both  (T/C  a .05  and  .02)  thickness/ chord  ratios 
(Figs.  12  to  14,  19).  The  velocity  components  agree  also  very  well  for  both  thickness/chord  ratios  (Figs. 

15,  16  and  21,  22).  From  drag  loop  comparisons  (Figs.  17  and  20a,  b)  at  near  mid  semi-span  only  slight 
differences  can  be  noted  in  the  peak  region.  From  figures  18  and  23  it  can  be  seen  that  there  is  also  good 
agreement  along  the  span  between  the  circulation  distributions. 

STRAKKD  WINO  Cases 

Only  results  for  the  T/C  - .02  case  are  available  for  datum  result  comparison.  The  agreement  between 
the  chordwise  pressure  distributions  on  the  inboard  wing  (Figs.  24,  25)  is  definitely  worse  than  for  the 
corresponding  thickness/chord  ratio  for  the  RAE  WINO  case.  On  the  outboard  wing  next  to  the  kink  the  agree- 
ment is  of  the  same  order  (Fig.  26).  Velocity  components  are  compared  for  the  two  stations  on  either  side 
of  the  kink.  The  agreement  is  fair  for  the  chordwise  velocity  component  (Figs.  27  and  29),  The  differences 
in  the  spanwise  velocity  component,  in  particular  at  the  inboard  section,  are  significant  (Figs  28,  30). 
Consequently  the  flow  directions  as  predicted  by  one  or  both  methods  are  in  error.  This  would  obviously  be 
of  some  importance  for  boundary  layer  calculations.  It  should  here  be  mentioned  that  there  is  some  indica- 
tion mat  Roberts*  results  are  not  fully  converged.  In  chapter  6.3.2  this  will  be  discussed  in  some  detail. 
Spanwise  circulation  distributions  again  show  good  agreement  (Fig,  31 )• 

NACELLE  Cases 

Datiim  results  for  the  nacelle  cases,  are  presented  in  figures  32  and  33.  For  these  eases  results  are 
also  a,railable  from  calculations  by  J.L.  Hess  (McDonnell  - Douglas),  These  results  were  obtained  using  a 
special  "higher-order"  axi -symmetric  program,  using  240  panels  chordwise.  From  figures  32,  33  it  can  be 
seen  that  there  is  a fair  agreement  between  all  three  methods.  Since  in  Hess*  method  no  circumferential 
discretisation  is  employed,  and  this  discretisation  contributes  sign  fi cantly  to  bo  the  error  in  internal 
flow  problems  (Hess,  Ref.  9),  this  solution  in  probably  the  best  of  the  datum  solutions  presented. 


5 .2  Comparison  of  "first-order"  Method  Hoaultn 

5.2.1  Introduction 

In  this  chapter  results  from  the  limit  - Semple  "sheets"  method  and  the  NLR-method  are  compared  graph- 
ically in  terms  of  chordwise  pressure  distributions,  drag  loops,  velocity  components  and  integrated  quant i- 
tieo  such  as  circulation  and  sectional  load*  It  may  be  noted  that  the  NLA-met had  in  considered  as  represen- 
tative of  the  first-generation  panel  methods* 

At  the  end  of  this  section  a few  NLH  results  are  also  compared  with  results  obtained  from  calculations 
carried  out  with  the  Hunt  - Semple  "lines"  method*  As  outlined  in  chapters  A.  2 and  A.A.l  these  methods 
employ  essentially  the  same  numerical  scheme. 

The  test  cases  were  run  with  prescribed  panel  distributions  defined  in  tables  1A  and  These  distri- 
butions are  considered  as  typical  engineering  distributions  for  most  "first-order"  panel  method  applica- 
tions in  which  detailed  pressure  distributions  are  required*  However,  it  should  be  noted  that,  although  60 
chordwise  panels  have  been  used  in  both  the  NLH  and  the  Hunt  —Semple  "sheets"  method  in  the  calculations 
discussed  in  this  section,  AO  chordwise  panels  are  usually  adequate  for  use  with  the  Hunt  - Semple  "sheets" 
method,  as  is  illustrated  by  the  results  discussed  in  chapter  6.2* 

5.2.2  Discussion 
RAF!  NINO  Cases 

The  chordwise  pressure  distributions  predicted  by  the  various  methods  for  the  RAE  W1N0  with  T/c  .IS 
and  a ('  were  found  to  be  practically  identical*  An  example  is  shown  in  filin'  A4.  It  can  be  seen  that 
for  this  non— lifting  case  agreement  between  the  fi  rst— order  NLR  results  and  Roberts*  solution  is  very  good. 

Results  for  the  RAJ*!  NINO  with  t/c  *15  and  a L°  are  presented  in  figures  i5  — A8.  The  pictuiv  n 

generally  similar  to  that  at  a 0 * However,  the  drag  loop  comparison  presented  in  figure  Ac  shows  some 

discrepancies  in  the  leading  edge  region  for  the  limit  - Semple  "sheets"  results.  Since  the  output  point; 
used  by  both  NLH  and  limit  - Semple  program  are  nominally  at  flat  panel  centroids,  the  difference  in 
values  near  the  loading  edge  must  be  dissociated  with  the  fact  that  NLR  effectively  optimise  their  chordwise 
amber  surface  doublet  variation  with  respect  to  sectional  local  lift*  Which  means  that  the  leading  edge 
sue t ion  pe.'ik  values  are  probably  overestimated  to  compensate  for  the  inescapable  At'  underesti mates  near 
the  trailing  edge,  associated  with  the  "lines"  modelling.  Thus,  the  accurate  NLR  drag  loop  peak  values 

would  appear  to  be  the  result  of  a rather  fortuitous  cancellation  of  independent  errors  associated  with  the 

/(’  and  the  0 values.  There  is  also  a noticeable  di  fferonce  in  the  spanwise  distribution  of  the  sectional 
load  (Fig.  A8,  note  that  a datum  solution  for  this  case  is  not  available).  The  reason  for  this  difference 
was  not  well  understood*  However,  the  span  wise  distributions  of  circulation  (Pig.  AY)  arc  virtually  iden- 
tical • 

For  T/h  .05  results  are  presented  in  figures  Ad  - ,)6*  For  the  chordwise  pressure  distributions 
there  is  again  good  agreement  between  the  various  solutions,  although  some  discrepancy  is  noticeable  in  the 
NLR  results  near  the  tip  (Fig*  41)*  The  Hunt  - Semple  results  again  differ  significantly  from  the  datum 
solution  in  the  leading  edge  region,  as  oan  be  seen  from  the  drag  loop  comparisons  presented  in  figure  42. 
Regarding  the  spanwise  velocity  component,  it  can  be  seen  from  figure  44  that  the  values  calculated  by  the 
NLR'  panel  method  are  in  error  over  a large  part  of  the  wing*  The  explanation  of  this  discrepancy  has  al- 
ready been  given  in  section  A. 3*2*  The  agreement  between  spanwise  distributions  of  circulation  and  load  is 
good  for  both  the  "first-order"  methods  (Fig.  45  and  46). 

Results  for  the  extremely  thin  RAE  WTNCJ  case,  i.o.  T/C  .0;'  are  given  in  figures  47  ~ 52.  Regarding 
the  chordwise  pressure  distributions,  it  can  be  seen  from  figure  47  that  there  is  poor  agreement  between 
the  NLR  results  and  the  datum  solution.  This  can  also  be  seen  from  the  drag  loop  comparisons  presented  in 
figures  48a,  b.  Details  of  the  pressure  distribution  in  the  leading  edge  region  are  again  poorly  predicted 
by  the  Hunt  - Semple  method*  It  can  be  seen  that,  in  the  NLR  results  for  this  cose  on  error  in  is  appa- 
rent, whereas  the  error  in  V has  increased  significantly  in  comparison  with  the  situation  for  T /C  .0*' 
(Figs*  4‘lf  50).  Also,  there  is  a significant  error  in  the  spanwise  load  distribution  predicted  by  the  NLR' 
method  (Fig*  52). 

From  the  observations  mod©  above  it  can  be  concluded  that  a "lines"  model  (as  employed  in  the  NLR 
panel  method)  generally  is  not  adequate  to  predict  aerodynamic  quantities  accurately  and  in  sufficient 
detai l for  wings  with  a thickness/ chord  ratio  below  .OS,  at  least  net  with  a panel  distribution  as  used 
for  these  calculations*  For  this  low  thiekneso/chord  ratio  only  the  circulation  can  be  calculated  with 
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reasonable  accuracy  using  ouch  a method. 

The  Hunt  - Semple  results  generally  agree  quite  well  with  Roberta  * datum  uolutionof  except  in  the  nose 
region.  This  is  possibly  explained  by  the  following.  Prom  figure  6cf  it  can  be  seen  that  the  use  of  the 
"optimiser"  reduces  the  source  strength  gradient  over  the  greater  part  of  the  chord,  but  not  in  the  imme- 
diate vicinity  of  the  loading  edge.  Consequently  the  discretisation  errors,  which  are  proportional  to  the 
source  strength  gradient,  can  still  significant  in  this  region. 

STKAKKD  WING  Cases 

Results  for  the  STRAKBD  WING  case  with  T/C  - .05  are  presented  in  figures  *>3  “ 64.  Prom  the  compari- 
sons of  pressure  distributions  at  three  inboard  and  three  outboard  stations  it  can  be  seen  that  the  agree- 

ment with  Roberts*  datum  solution  is  reasonable,  both  inboard  and  outboard  for  both  methods  (Pigs.  53  — 5®), 
the  worst  errors  being  shown  by  the  NLR  results  just  inboard  of  the  crank  station  (Pig.  55)*  Spanwise  and 

chordwise  velocity  components  are  presented  for  the  stations  at  either  side  of  the  kink  in  figures  59  -62. 

It  can  be  seen  that  the  error  in  at  the  inboard  section  is  significant  for  the  NLR  results  (Pig.  60).  In 

contrast,  the  Hunt  - Semple  results  exhibit  a much  smaller  error.  This  error  is  probably  due  to  the  fact 
that  in  the  Hunt  - Semple  method  a 2-D  "optimiser"  is  used  for  the  determination  of  the  internal  doublet 
distribution  (see  3*3*l)»  whereas  the  flow  at  this  station  is  highly  3-D. 

Por  T/c  - .02,  results  are  presented  in  figures  65  “ 73*  Inboard,  the  pressures  are  again  predicted 
reasonably  well  by  both  methods f however,  the  NLR  errors  are  larger  than  those  of  the  Hunt  - Semple  method. 
At  the  outboard  station  next  to  the  kink,  the  errors  in  the  NLR  results  are  similar  to  those  in  the  mid- 
semispan  results  for  the  corresponding  RAE  WING  case  (Pig.  47 # T/C  - .02).  Regarding  the  comparisons, 
it  can  be  noted  from  figure  68  that  at  the  inboard  station  the  error  in  has  increased  slightly  for  both 
methods  in  comparison  with  the  corresponding  plot  for  T/C  - ,05  (Pig*  59)* 

At  the  outboard  section  next  to  the  kink  the  situation  is  similar  to  the  RAP'  WING  case  with  T/c  « .02  i.e. 
the  Hunt  - Semple  results  agree  well  with  the  datum  solution,  whereas  the  NLR  results  are  seriously  in  error. 
Regarding  the  comparisons  the  situation  is  quite  different.  Inboard  both  methods  are  seriously  in  error 
(Pig.  69)  j n.b.  two  datum  solutions  are  shown.  However,  figure  69  indicates  some  uncertainty  in  datum  solu- 
tion accuracy.  At  the  outboard  section,  however,  the  results  again  agree  reasonably  well  with  Roberta*  da- 
tum solution. 

The  sectional  load  distributions  are  presented  in  figure  73*  Apart  from  small  differences  on  the  inboard 
part  of  the  wing,  the  Hunt  - Semple  results  agree  very  well  with  Rubbert’s  solution.  Por  the  NLR  results 
the  situation  is  different,  in  that  on  the  outboard  wing  the  comparison  is  similar  to  that  for  the  RAE  WING 
with  T/C  « .02,  i.e.  the  error  in  the  load  is  significant,  whereas  on  the  inboard  part  of  the  wing  the 
agreement  with  Rubbert*s  results  is  relatively  good . However,  the  distinct  jump  in  the  NLR  results  across  the 
crank  seems  unrealistic. 

Summarizing  the  "first-order"  method  comparisons  for  the  STKAKKD  WING  eases  it  may  be  said  that,  for 
the  results  on  the  outboard  wing  the  same  cone 1 us ions  apply  as  those  al ready  given  for  the  RAE  WING  cases. 
Regarding  the  inboard  part  of  the  wing,  errors  in  the  predicted  spanwise  velocity  component  become  very 
apparent,  particularly  for  the  T /C  ,0.'  case.  Por  the  Hunt  - Semple  results,  this  may  possibly  be  due  to  the 
2-D  nature  of  the  **optimi ser**  used  in  this  method  (set  chapter  '.'.l).  However,  it  should  be  mentioned  that, 
apart  from  "sheets**  modelling  being  inherently  better  than  that  of  the  "lines"  type*,  the  Hunt  - Semple  meth- 
od does  allow  the  user  to  modify  the  results  of  th.  usual  2— D Optimiser",  if  it  is  felt  advantageous  to  do 


NACELLE  Cases 

Chordwise  pressure  distributions  for  the  NACp.LLR  case  with  chord/exitdiamcter  ratio  1.0  and  an  angle 
of  attack  of  5°  are  compared  at  three  circumferential  stations,  and  are  shown  in  figures  74  “ 76.  At  all 
three  stations  the  Hunt  - Semple  results  agree  well  with  Rubbert*s  datum  solution  on  the  outside  of  the 
nacelle.  On  the  innorside,  some  discrepancies  can  be  seen.  Generally,  the  NLR  results  do  not  show  such  good 
agreement  as  those  from  the  Hunt  - Semple  method. 

In  figure  77  results  are  compared  for  chord/exi t diameter  ratio  3.333*  at  an  angle  of  0°.  The  flow  for 
this  type  of  configuration  can  be  characterized  as  being  more  of  the  ’•internal**  type  than  in  the  case  dis- 
cussed above.  It  can  be  seen  from  figure  77  that  inside  the  nacelle,  both  "first-order"  methods  arc 
seriously  in  error.  These  errors  can  bo  associated  with  the  fact  that,  for  internal  flows  of  this  type,  the 
conservation  of  mass  is  not  satisfied,  the  leakage  being  due  to  "fi rot-order"  modelling  deficiencies. 


3*2*3  Comparison  of  NLR  Results  with  Hunt  - Semple  "lines"  Method  Results 

A comparison  is  presented  to  illustrate  the  type  of  differences  that  may  occur  when  calculations  are 
carried  out  with  two  methods,  that  basically  employ  an  identical  numerical  model  with  regard  to  surface 
approximation  and  the  use  of  a discrete  vortex  system  (i*c*  a "lines"  model). 

Differences  between  the  results  than  are  only  due  to  a different  chordwise  variation  of  the  internal  vortex 
strengths  (note  the  use  of  the  "optimise!*"  in  the  Hunt  - Semple  "lines"  method,  see  chapter  3*3*1 ) and  the 
formulation  of  the  Kutta-condition  ("calculated"  in  the  Hunt  - Semple  method). 

Results  for  the  RAE  WING  with  T/C  ~ ,03  are  given  in  figures  J8  -82,  Prom  figure  78  it  can  be  seen 
that  the  pressure  distributions  at  the  section  near  mid— semispan  agree  to  about  the  same  level  of  accuracy 
with  Roberts’  results.  This  applies  also  to  the  comparisons  of  the  velocity  distributions  (Pigs,  /9,8o) • 
Agreement  between  the  circulation  distributions  is  again  very  good (Pig*  8l),  Prom  figure  82  it  can  be  noted 
that  the  Hunt  - Semple  "lines”  method  tends  to  underpredict  the  load  slightly. 

6.  SWSlTIVm  OP  THE  RESUUFS  TO  THE  NUMBER  OP  PANELS  USED 
6.1  Introduction 

In  this  section  some  results  are  presented  and  discussed  which  were  obtained  from  calculations  carried 
out  using  different  numbers  of  panels.  For  both  the  NLR  panel  method  and  the  Hunt  - Semple  "sheets”  method 
a common  panel  distribution  was  prescribed  in  which  only  30  panels  chordwise  were  used.  Spanwise  the  dis- 
tribution is  the  same  as  for  the  preceding  comparisons  (see  table  IB  and  2).  Also  there  are  some  results 
available  from  the  methods  of  Rubber t and  Roberts’  which  were  obtained  for  fewer  panels  than  were  used  to 
obtain  datum  accuracy.  For  both  these  methods  the  panelling  used  was  left  entirely  to  the  choice  of  the 
user. 


6,.'  Results  from  "first-order"  Methods 
RAE  WING  Cases 

Results  from  the  NLR  method  for  the  RAE  WING,  with  T/C  .13  are  presented  in  figure  83  — 8‘«.  These 
indicate  that,  for  a relatively  thick  wing,  30  chordwise  panels  are  sufficient  to  give  good  accuracy. 

Some  results  for  the  RAE  WING  with  T/c  .03  are  given  in  figures  86  - 89.  Prom  the  pressure  distri- 
butions near  mid-semispan  it  appears  that,  in  order  to  get  a good  quantitive  prediction  of  the  pressure;-., 
30  chordwise  panels  suffice  when  the  Hunt  - Semple  "sheets"  method  is  used.  For  a "lines"  method  like  that 
of  NLR,  30  chordwise  panels  are  not  sufficient  to  predict  the  pressure  distribution  in  detail.  The  span- 
wise  velocity  component  at  mid-semispan  is  given  in  figure  87*  Apparently  the  results  obtained  from  1 1. « • 
Hunt  - Semple  "sheets"  method  with  30  chordwise  panels  are  practically  as  accurate  as  the  (60  x 12)  re- 
sults. The  NLR  results  show  a worsening  of  the  type  of  error  already  present  in  the  (60  x 12)  results,  as 
the  number  of  panels  is  decreased.  However,  for  the  circulation  distribution  (Fig*  88),  the  agreement  with 
Roberts*  datum  solution  for  both  the  Hunt  - Semple  "sheets"  method  and  the  NLR  method  is  seen  to  be  re- 
markably good  for  only  30  chordwise  panels.  Figure  89  illustrates  that,  in  order  to  predict  the  spanwise 
loading  with  acceptable  accuracy  the  NLR  method  must  be  used  with  many  more  than  30  chordwise  panels. 

Results  for  the  T/C  .02  case  arc  given  in  figures  90  - 93*  Prom  figure  90  it  can  be  seen  that  the 
pressure  distribution  calculated  with  the  Hunt  — Semple  "sheets"  method,  and  employing  only  30  chordwise 
panels,  still  show  a good  accuracy.  For  this  case,  results  from  the  NLR-method  are  available  with  up  to 
90  chordwise  panels.  The  NLR  results  using  a (90  x 12)  panel  distribution  show  a distinct  improvement  in 
comparison  with  the  corresponding  (60  x 12)  results.  For  this  case,  the  (30  x 12)  NLR  results  are  very  in- 
accurate. The  spanwise  velocity  component  comparison  is  presented  in  figure  91,  The  Hunt  - Semple 
"sheets"  results  are  accurate  for  both  the  panel  distributions  used.  On  the  other  hand  the  NLR  results 
show  that,  increasing  the  number  of  chordwise  panels  up  to  90  does  not  lead  to  a significant  improvement 
iri  accuracy  of  the  spanwise  velocity  component.  However,  from  figure  92,  it  can  be  seen  that  even  for  such 
a low  thickness/ chord  ratio  (i.e.  .02)  the  circulation  is  still  predicted  wi th  good  accuracy  by  all  meth- 
ods. With  regard  to  the  sectional  load  distributions  presented  in  figure  93,  the  situation  is  quite  diffe- 
rent. The  Hunt  - Semple  "sheets"  method  is  capable  of  predicting  the  load  distribution  with  good  accuracy^ 
even  using  only  (30  x 12)  panels.  However,  the  NLR  results  again  show,  in  particular  for  the  (30  x 12)  dis- 
tribution, a rather  poor  agreement  with  the  datum  solution.  Furthermore,  it  can  be  noted  that  even  the  use 
of  the  ( K)  x 12)  panel  distribution  is  not  sufficient  to  obtain  acceptable  accuracy. 


STRAKKD  tfINCI  Case s 

In  Mgur«i  94  - 100,  results  are  presented  for  the  two  STRIKED  W INQ  i'«h«h  i.e.  for  T/C  - ,0*>  and  .02. 

For  T/C  - ,0*>,  figure  94  shows  the  ohordwise  preaaure  distributions  plotted  at  the  inboard  section 
ne*t  to  the  orank|  thie  section  le  considered  to  be  the  most  difficult  section  for  thie  particular  confi- 
guration. Tile  agreement  of  the  ( AO  x 12)  results  ran  gee  from  good  (Hunt  - temple  "sheets"  method)  to  rea- 
sonable (NLH-method) , Tlie  good  agreement  between  the  Nl.K  ( 10  x 12)  reaultu  and  the  datum  eolution  le  un- 
expected in  view  of  the  relatively  poor  agreement  for  the  RAE  NINll  case  with  T /C  - .OS  (eee  Fig.  86). 

The  apanwiee  velocity  components  for  thie  oaee  are  compared  in  figure  9S,  It  >-aii  be  aeon  that,  though  the 
preaaure  dletributione  agree  fairly  well  at  thie  section,  the  reeulta  do  notf  the  Hunt  - Semple  "sheets" 
reeulte,  though  in  error,  being  considerably  better  than  the  beet  from  the  NLtt-m  stUnt.  The  Her  t tonal  load 
diet  nbutione  obtained  from  both  the  NLR  and  the  Hunt  - Semple  met  hod  e are  presented  in  figure  >/.  It  ran 
be  seen  that  a decrease  of  the  number  of  panela  need,  to  10,  oorreaponda  to  a uigni  fi  rant  derreaae  In  the 
accuracy  of  the  calculated  load  diatrlbution  for  the  NLR-method,  whereas  the  agreement  between  the  Hunt  - 
Semple  reaultu  for  ( 10  x 12)  and  (60  x 12)  reunite  is  very  good  (n.b,  no  datum  aolutlor.  for  thie  cane  le 
available).  In  figured  k>8  - 100  reaultu  are  preuented  for  the  extremely  thin  STRAKKb  WING  caae,  i.e. 

T/C  - .02*  The  ohordwise  preaaure  di  at  ri  but  i on  a for  the  moat  critical  auction  i.e.  the  auction  juet  inboard 
of  the  crank  are  given  in  figure  ‘>8.  Tleaue  note  that  two  extra  aeta  of  results  have  been  int  i\iduoed,  which 
were  obtained  uaing  the  NLR  panel  method.  It  can  be  seen  that  for  the  NUf-method  10  chorvlwiae  i>anela  are 
definitely  inadequate  to  predict  the  preaaure  diatrlbution  with  auffioient  detail.  Further,  it  can  be  no- 
ted that  an  increase  from  60  to  ohordwise  panela  hardly  changed  the  preaaure  diatrlbution.  An  increaae, 
how  ever,  of  the  number  of  apanwiae  strips  fivm  12  to  18  (n.b.  6 extra  at  rips  in  the  vicinity  of  the  crank 
were  uaed)  1 cade  to  some  improvement . 

A span  wise  velocity  component  comparison  u given  m figure  09.  Here  it  can  also  be  aeen  that  an  increaae  of 
the  number  of  chorvlwiae  panela  hardly  improved  the  level  of  accuracy.  In  fact,  increasing  the  number  of 
apanwiae  utripu  clearly  given  the  moat  algnificant  improvement  in  accuracy,  Sectional  load  distributions 
are  cornered  Hi  figure  100.  The  agreement  of  the  Nl.R  reeulta,  obtained  using  10  chorvlwiae  panela,  with 
Nubbert9*  datum  aolution  la  rather  poor.  Further,  it  can  be  aeen  that  even  when  uaing  a (k'  \ 12)  distri- 
bution the  predicted  load  ia  a till  about  l 't  in  error. 

o.i  Results  from  ••higher- order"  method* 

6.1.1  Rubbert’a  method 

From  Rubbert*s  method  reaultu  are  available  for  the  RAK  NINO  case  with  T,/f  - ,0S,  obtained  using  two 
"engineering"  panel  vtiet ributiona,  i.e,  one  with  22  ohordwise*  panels,  and  12  apanwiae  strips  and  one  in 
which  only  12  chorvlwiae  panels  are  employed,  but  again  with  12  apanwiae  strips.  In  practice,  a diatrlbution 
with  22  chord wise  panels  and  only  8 apanwiae  strips  could  be  used  to  give  reasonably  accurate  results  for 
such  a simple  wing  configuration • 

In  figure  101  the  ohordwice  pressure  di at ributiona  near  mid-semi  a pan  are  presented  and  compared  with 
Roberts*  datum  solution.  It  can  be  seen  that  the  (22  x 12)  distribution  gives  very  accurate  results.  Using 
a (12  x i ) distribution,  which  is  characterized  by  a coarser  panelling  near  the  leading  edge,  loads  to  a 
loss  of  accuracy  in  this  region.  From  figure  102  it  can  be  seen  that  the  apanwiae  velocity  component  is 
predicted  very  well  even  when  only  22  ohordwise  panels  are  used.  The  sectional  load  dial ributiona  plotted 
in  figure  101  show  that,  although  the  chord wise  pressure  distributions  are  not  correct  in  detail,  the  sec- 
tional load  is  still  predicted  with  good  accuracy!  even  when  the  (12  x 12)  distribution  is  used. 

6.1.2  Roberts *-method 

The  sensitivity  of  the  results  from  Robert s '-method  to  the  number  of  ivuiela  used  is  demonstrated  for 
both  the  "thin"  RAK  NINO  and  the  "thin”  ST  RAK  ED  NINO  case.  For  the  RAK  NINll,  results  are  available  from 
calculations  carried  out  with  two  different  panel  distri buttons,  i.e.  one  with  27  ohordwise  panels,  and 
> spanwise  strips  and  a second  one  which  has  l4)  ohordwise  panels  and  only  6 apanwiae  strips. 

In  figure  104  the  ohordwise  pressure  distributions  near  mid— semi  span  are  presented  for  the  RAK  WING  case. 
The  agreement  between  the  datum  solution  and  the  results  from  both  panel  distributions  is  very  good  over 
a large  part  of  the  chord.  From  figure  104  it  can  be  noted,  however,  that  when  the  number  of  panels  used 
is  decreased,  the  pressures  show  a rapid  increase  of  error.  It  could  even  be  questioned  whether  the  datum 
solution  ia  "fully  converged"  near  the  trailing  edge,  Spanwise  velocity  components  are  presented  in 
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f l gill**  10*'*  dolirPAlly,  the  Agreement  between  the*  0O  l U t i ono  10  Very  good,  A l t hough  the*  name  teildoliov,  1 ,«- , 
the'  rapid  lUOVvAOe  of  iH'lVr  Hear  the’  trailing  edge  e .in  he*  Ooell  , 0 1 I'Oul  At  | on  ill «» t 1*1  t»Ut  1 .’III  Art*  presented  111 
flgUfe  I0f-,  A p I'M  1**11 1 ly,  t tie*  Utt*  Of  fewer  pAlielo  reoultn  ill  Alt  O VePpivd  1 0 t l Oil  of  t tie’  » l IV  u 1 <4 1 1 Oil  along  ttle* 
OpAll  * 

For  the  1* t\AK Kl'  WINd  moulto  Are*  available  for  the*  tA*  .0  OAUe*  Choi'dWl  tie*  ptvndUPe  ill  jit  II  tuit  1 oil.: 

-»t  the-  Ulb.'Alxl  O t A t l OH  Holt  to  ttle*  O rank  A IV  p IVO  e’ll  t ed  ill  fl^UV  10/,  A deotvatie  (ill  ttle*  \*  l . * 1 1 1 1 t y O f tile 
• Alik  ) of  ttle*  III  UUt'e»  I*  Of’  Ul'AllWl  tie*  tit  l*i  pel  laird,  1 rate  to  fe*  l A t 1 Ve*  t y UIIIaII  o hail  gen  ill  ttle*  pleO.O.UIV  dlntrlbu”* 
tl.’li  At  ttll.  aootloil,  A fl  I tfll  1 A I*  > ompAPi  *»*VI  in  ob  t A-l  Med  for  ttle*  tl|  All  W)  O t*  Velooity  . *«'lilpOfl  t*ll  t f (Fig.  1 1>M  ) , 

Ttle*  Oi  I 11 1 At  I Oil  ill  il  t 1*1  t'U  t i oil  Apparently  in  tlAlMly  altered,  AO  >Ail  tv  neell  t‘ 1*0111  figure  U>), 


AOiltKAiT  VKICUItf  iVMIOTATWW  PI  MW*. 

Ill  thin  vtlAptfl*  A opiate*  Attempt  ill  made  to  t*  VA 1 UA  t e the  relative  Ille’Pltti  of  ttle*  Uiethodt  ooilipared,  III 
ten*.:  of*  ttle  output  At  l on  A l effort  required  to  ohtAlli  A oertalll  level  of’  AOOUrAOy,  Ail  men  t l oiled  al  P»  ady  ill 
tiApter  l,  tht:  pone.*  two  pro  Memo,  ’'iretly,  ttie  o a l on  l a t i ono  with  t tie  VAPioilti  programo  were  . Arriril  out  Oil 
different  e output  era*  Ttilri  niAjv  e o it  i iiipono  l M r to  ooitipAre*  ttle  reopeotive  oomputatloii  t i men  in  a itlivot  w Ay, 

• •«*  on.tlyv  t tie  level  of  ttie*  aoouraoy  muo  t he  defined  in  oome  appropriate  manner. 

I'he  firnt  problem  tiAo  been  approaohed  hy  oral  lug  aII  oomputmg  1 1 men  to  ttie)  time  refill  reset  on  a refer- 
er.  r eomputerj  for  thin  purpooe  ttie  OlX-  f’f’OO  tiAti  been  ohooeit . "tie  oeallng  far  torn  uoed  (nee*  table  4)  tiAVe 
been  eo  t \ n.it  ed  by  the  iiutlVivtUAl  pAPt  i ot  p.ui  t n in  t til  u i ii  veu  t 1 ga  t l on  , It  will  be  olear  that  thenr  fAotora 
dtould  be  l n t erptv  t ed  wltti  OAPe,  For  refereiioe  purpooeo , oome  of  ttie  AotiiAl  ooiliputatlon  timer  are  preoeuted 
in  table  , Kor  three  of  the  teot  oaneu  two  ill  fferent  timeo  are  gi  veil , Ktrntly,  ttie  time  for  ttie  . output  er 
on  whioh  ttie  oaoe  wao  AotiiAlly  nut  and,  t<e  tween  braoketo,  ttie  entimated  time  for  a MX'  f'C»00,  uomg  the  noal- 
ing  f.i.  toPii  fivm  table  4, 


I'he  level  of  a*'.*up.i>'v  h an  been  defined  by  me  at  in  of  a no— railed  Hl,  —error  norm",  Thin  l,  error  norm  V' 


t given  over  oome  interval  < by; 
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In  parti,  ul  ir  the  following  .h  o.  re  t i nod  fair  of  e.pi  (.‘0  ttao  been  lined; 
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I'he  c-riv't*  norm  K a.  •noiMliif'  to  e.pi  (*)1  ha.  been  determined  foi*  preooure  dint  n but  1 ono  (v^  s Al'  ) .and  op.uiwme 
veto,  i ty  •'ompotieii  t dl  o t rt  but  i v'no  (b  * AV^)f  wtierevn  lo  the  differein'e  between  Koberto*  datum  notation, 

alnulated  at  N ^ - KX'l  e.pt  idiot  .ui  t poutto  on  tin*  interval  < { * ohord),  and  any  other  nolution.  The  valueo 
it  the  poliito  i were  found  throu^ti  a oplute  i n t arpo  l a t ion  prooedure, 

:’oi«e  reoultn,  obtained  from  the  I'rooedure  outlined  above,  are  plotted  m figured  1 10,  III,  Krror  nonuo 
for  rhordwi  e preooure  dl ot n but i ono  uid  opanwine  velooity  dt o t rt but  tone  are  oonnidered  at  ml d—eemi opan 
(H  - .** 4 *0  for  two  h'AK  WlNil  oaoen  (T/0  - ,0*>  uid  and  at  4 - ,.‘lv>  (juot  inboard  next  to  the  orank'for 

the  two  dT  h’AK  Kb  WINd  o Aden  • 


Ke^nPxtiiig  the  h'AK  WlNvl  oaoeo,  it  oon  be  net'll  from  figure  110a,  o that  for  all  methodo  the  erivr  m v' 
ln.Mv.iueo  n i gni  f t .Ail  t ly  when  the  thtoknenn  ohord  ratio  deoreauen  from  ,0*'  to  ,0.',  Further,  it  ma^V  be  noted 
that,  for  a given  oomput  at  t ona  1 effort  the  error  m Rlibbert*a  reaulto  in  omaller  than  that  of  all  other 
methi'itn.  Flirt  her  | it  oaii  be  neon  that  he  dlfferenoe  in  error  norm  E between  the  NLR  and  the  Hunt  - Oemple 
**nheetow  reoultn  lo  aotnewhat  oiualler  than  migtit  be  expeoted  from  the  oorreopondmg  preooure  diotribution 
ptotn.  Thin  in  pooolbly  due  to  (tie  relatively  large  erroro  in  the  Hunt  - Sample  HoheetoM  reoultn  preoent 
in  the  leading  edge  region  • Aloo,  it  oaii  be  oeen  in  figure  llOo,  that  the  error  norm  for  Rubbert  *h  datum 
nolution  and  for  Roberto*  engineering  oaoe  (.'/  x v>l  p.anelo  in  of  the  name  order  of  magnitude  but  that  the 
oorreopondmg  oomput  at  ion  timeo  differ  algni  ft  rant  ly  • 

Regal'd  mg  the  opanwine  velooity  oomponento  the  oituation  in  very  similar* 

In  figure  III  reoultn  of  the  name  type  are  preoented  for  the  STKAKKl'  WINd  oaoeo  at  the  oeotlon  Juat 
inboard  of  the  orank*  denerally  the  name  oonoluoioue  apply  ao  have  already  been  drawn  for  the  KAK  WINd 
oaoeo.  In  addition,  It  oaii  be  noted  from  the  Nl.R  reoultn  that  an  moreaoe  of  the  number  of  opanwine  ntnpn 
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(with  the  number  of  chord  wise  panels  kept  fixed)  gives  a signi  ficantly  better  prediction  of  the  spaiiwise 

velooity  component  than  a chord wise  increase  with  the  spanwise  distribution  fixed. 

'Hie  C -error  nona  results  from  the  NLR-raothod  for  the  RAE  NINQ  cases  are  plotted  as  function  of  the 
thlokness/chord  ratio  in  figure  11.' • This  figure  shows  the  rapid  increase  of  the  errorf  when  only  30  x 1? 

panels  are  used,  for  decreasing  thickness/chord  ratio.  Also  it  can  be  seen  clearly  that  when  60  x 1.'  panels 

are  used,  good  results  can  be  obtained  employing  a "lines"  model  (as  used  in  the  NLIHnethod)  for  relatively 
thick  wings,  i.e.  with  the  thioknoss/chord  ratio  not  below  .0*). 

Finally,  it  is  worthwhile  bringing  attention  to  the  fact  that  both  the  NLH  and  the  Hunt  - Semple  pro- 
grams are  very  general  in  terms  of  application  scope  and  are  extremely  "user  oriented"  in  terms  of  input/ 
output  facilities  etc,.  Because  of  their  generality  of  scope  and  design  for  ease  of  use,  such  programs  are 
always  liable  to  be  less  efficient,  for  any  particular  simple  case,  than  programs  having  a smaller  practi- 
cal all-round  utility  value. 

8.  CONCLUSIONS 

A comparative  study  has  been  made  of  the  capabilities  and  efficiency  of  several  so-called  "panel  meth- 
ods", with  respect  to  the  prediction  of  aerodynamic  quantities  such  as  pressure  di strlbut i one,  velocity 
distributions,  circulation,  and  integrated  quantities  such  as  sectional  load.  The  methods  currently  compa- 
red are:  l)  Roberts*  Sp  1 1 n e-N eumann  Program,  developed  at  HAe  (Weybndge),  i i ) Hc«  ingV  Interim  Higher  Order 
Method,  developed  by  Rubbert  et , al.,  iii)  the  Hunt  - Semple  "sheets"  method  and  iv)  the  NLK  Panel  method* 

The  configurations  for  which  calculations  have  been  carried  out  are:  l)  n set  of  simple  swept  wings 
(KAK  NINO  cases),  ii)  a set  of  wings  with  strake  (BTRAKKP  NINO  cases)  uni  iii)  a set  of  ring-wings  or  na- 
celles ( N AC  F.LLV  c as es  ) . 

High  accuracy  (datum)  solutions  were  provided  by  the  methods  of  Roberts  and  Rubbert, 

Regarding  the  wing  test  cases,  the  following  has  beer  observed: 
l)  It  appears  that  generally  the  agreement  between  the  datum  solutions  from  Roberts*  "thi rd-order**  Spline- 
N eumann  program  .uni  Rubbert  *s  "pseudo  second-order"  method  is  very  good.  However,  on  the  inboard  part  of 
the  STRAKKP  NINO  ease,  with  thi oknoss/ohord  ratio  .0*',  some  dlscrepanci es  can  be  noted  in  particular  for 
the  spanwise  velocity  component.  Unfortunately  it  cannot  be  decided  which  of  tin'  solutions,  if  any,  i;  the 
"correct  one", 

ii)  It  appears  that  first  generation  panel  methods  employing  discretized  vortex  sheets  (like  the  NI.K  panel 
method,  the  HUB  method,  the  Boeing  TKA  .'10  program  and  the  Hunt  - Semple  program  using  the  "lines"  option’ 
are  adequate*  for  predicting  aerodynamic  quantities  for  relatively  thick  wings  (i.e,  thickness/ chord  ratio 
not  below  .0*'). 

iii)  The  limit  - Semple  program  using  the  "sheets"  option  offers  an  improved  capability  for  "thin  wing" 
applications!  and  is  also  relatively  more  efficient  for  calculating  "thick  wings", 

iv)  Considering  relative  merits  in  terns  of  accuracy  versus  associated  computation  time,  with  different 
numbers  of  panels  being  used  for  each  method,  strongly  suggests  that  the  Boeing  Interim  Higher  Order  Method 
is  the  most  efficient.  However,  it  should  be  noted  that  some  of  the  methods  have  evolved  considerably  since 
the  currently  compared  results  were  actually  calculated  (i.e.  in  l>J/’t'). 

Nith  regard  to  the  nacelle  test  cases  the  following  has  observed: 
i)  Discrepancies  can  be  noted  between  the  datum  solutions  for  all  the  nacelle  cases  considered,  in  particu- 
lar on  the  innerside* 

ii  ) On  the  innerside  of  the  slender  NAOKM.K  case  (chord /exit  diameter  ratio  l,  III),  there  is  roor  agreement 
with  tiie  datum  solution  for  both  "fi rot— order"  methods  (i.e.  NI.K  and  Hunt  - Semple  "sheets"  method)*  How- 
ever, it  may  be  inferred  that  these  discrepancies  are  not  so  much  due  to  "lines"  versus  "sheets"  modelling 
differences,  but  rattier  to  the  chosen  ohordwise  shapes  of  doutiUetty  variation  on  the  mean  surface  used  in 
conjunction  with  constant  strength  source  panels  on  the  geometric  surface. 
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TABLE  1A 

Chordwiae  definition  of  fixed  panel  distribution  (60  panels) 


NACA  POUR-DIGIT  SERI  BE 

t/c 

- 0.020000 

X/C 

z/c 

0 

0.0 

0.0 

1 

0.002084 

0.001329 

2 

0.008319 

0.002601 

3 

0.018656 

0.003808 

4 

0.033014 

0.004941 

5 

0.051281 

0 ,00S 484 

6 

0.073311 

0 .006937 

7 

0.098933 

0.007774 

8 

0.127444 

0.008489 

9 

0. 160117 

0.004071 

10 

0.145149 

0.004616 

11 

0.232413 

0.009816 

12 

0.272463 

0 .009975 

13 

0.315035 

0.009994 

14 

0.358797 

0.009882 

15 

0 .403405 

0.009647 

16 

0.450004 

0.009301 

17 

0.496733 

O.OO8858 

18 

0.543724 

0.008330 

14 

0.600606 

0.007733 

20 

0.637012 

0.007080 

21 

0.682576 

0.006385 

22 

0.726441 

0 .005661 

23 

O.764757 

0.004920 

24 

0.810687 

0.00417 3 

25 

0.849411 

0.003432 

26 

0.885622 

0.002709 

27 

0.914038 

0.002013 

28 

0.944394 

0.001358 

29 

0.476452 

0.000753 

30 

1 .000000 

0.000210 

I NACA  FOUR-DIGIT  SERIES  1 

t/c 

. 0.050000 

x/c 

z/c 

0 

0.0 

0.0 

1 

0 .002084 

0.003322 

2 

0.008319 

0 .006502 

3 

O.OI8656 

0.004520 

4 

0.033014 

0.012353 

5 

0.051281 

0.014471 

6 

0.073311 

0.017343 

7 

0.098933 

0.019436 

8 

0.127944 

0.021223 

9 

0.160117 

0 .022679 

10 

O.I95199 

0.023787 

11 

0.232913 

0.024540 

12 

0.272963 

0.024436 

13 

0 .315035 

0.024986 

14 

0.358797 

0.024705 

15 

0.403905 

0.024118 

16 

0.450004 

0.023253 

17 

0.496733 

0.022144 

18 

0.543724 

0.020825 

19 

0.5 90606 

0.019333 

20 

0.637012 

0 .017701 

21 

0. 682676 

0.015964 

22 

0.726941 

0.014153 

2? 

0.769757 

0.012300 

24 

O.8IO687 

0.010433 

25 

0.849411 

0.008581 

26 

0 .885622 

0.006772 

27 

0.919038 

0.005033 

28 

0.949394 

0.003394 

29 

0.976452 

0.001882 

30 

1 .000000 

0 .0005 25 

NACA  FOUR-DIGIT  SERIES 

t/c 

- 0.150000 

x/c 

z/c 

0 

0.0 

0.0 

1 

0.002084 

0.009967 

2 

0.008319 

0.019506 

3 

0.018656 

0.028561 

4 

0.033014 

0.037060 

5 

0.061281 

0.044914 

6 

0.073311 

0 .052028 

7 

0.098933 

0.058308 

8 

0.127944 

0.063668 

9 

0.160117 

0 .068036 

10 

0.195199 

0.071362 

11 

0.232913 

0.073620 

12 

0.272963 

0.074809 

13 

0.315035 

0.074957 

14 

0.358797 

0.074115 

15 

0.403905 

0.072353 

16 

0 .450004 

0.069759 

17 

0.496733 

0.066432 

18 

0.543724 

0.062476 

19 

0 .590606 

0.057998 

20 

0.637012 

0.053102 

21 

0.682576 

0.047891 

22 

0.726941 

0.042460 

23 

0.769757 

0.036900 

24 

O.8IO687 

0.031300 

25 

0.849411 

0.025744 

26 

O.885622 

0.020316 

27 

0.919088 

0.015100 

28 

0.944394 

0 .010182 

29 

0.976452 

0.005646 

30 

1 .000000 

0.001576 

TABLE  IB 

Chordwiee  definition  of  fix?d  panel  distributions  (30  panels) 


NACA  FOUR-DIGIT  SERI  S3 

t/c 

* 0.050000 

x/c 

z/c 

0 

0.0 

0.0 

1 

0.006234 

O.OO566I 

2 

0.025317 

0.010957 

3 

0.057991 

0.015766 

4 

0.105167 

0.019865 

5 

0.167863 

0.022962 

6 

0.246917 

0.024722 

7 

O.3422oq 

0.024847 

8 

0.451964 

0.023211 

9 

0.570710 

0.019986 

10 

0.690027 

0.015668 

11 

0.799534 

0.010951 

12 

0.889014 

0.006599 

13 

0.950584 

0.003329 

14 

0.984054 

0.001448 

15 

1.000000 

0.000525 

NACA  FOUR-DIGIT  SERIES 

t/c 

■ 0.150000 

x/c 

z/c 

0 

0.0 

0.0 

1 

0.006234 

0.016982 

2 

0.025317 

0 .032872 

3 

0.057991 

0.047297 

4 

0.105167 

0.059596 

5 

0.167863 

0.068887 

6 

0.246917 

0.074165 

7 

0.342298 

0.074541 

8 

0.451964 

0.069633 

9 

0.570710 

0.059958 

10 

0.690027 

0.047003 

11 

0.799534 

0.032853 

12 

0.889014 

0.019796 

13 

0.950584 

0.009986 

14 

0.484054 

0.004344 

15 

1.000000 

0.001575 

NACA  FOUR-DIGIT  SERI  IE 

t/c 

. 0.020000 

x/c 

z/c 

0 

0.0 

0.0 

1 

0.006234 

0.002264 

2 

0.025317 

0.004383 

3 

0.057991 

0.006306 

4 

0.105167 

0.007946 

5 

0.167863 

0.009185 

6 

0.246917 

0.009889 

7 

0.3422  >8 

0.009939 

8 

0.451964 

0.009284 

9 

0.570710 

0.007994 

10 

0.690027 

0.006267 

u 

0.7  ,,r'34 

0.004380 

12 

0.889014 

0.002639 

13 

O.950584 

0.001331 

14 

O.984OS4 

0.000579 

15 

1.000000 

0.000210 

Airfoil  section:  Symmetrical  NACA  TOUR-DIGIT  AIRFOIL 

Z/C  . ^ (.. 'l6'K>(x/c)2  - .12600( X/C ) - .35l60(x/c)2  + ,28430(X/C)3  - . 10150(x/c)4 


KAK  WINO,  a plan  form 

n Valued  or 
Panel  Edit'd 

n Valued  of  Output 
St at 1 one  ( Cen t roi do ) 

0.0 

0.0043 

0.049 

0.0788 

o.ior 

0.1466 

0.186 

0.2314 

0.079 

0.3307 

0.184 

0.4380 

0.494 

0.5489 

0.606 

0.6683 

0.713 

0.7604 

0 .810 

0.86 Or 

0.893 

0.9044 

0.967 

O.978O 

1.0 

STKAKK1'  W1N0  plan  form 


9 Values  of  n Values  of 
Panel  Kill's  Output  Stations 
(Cent  rcids) 


0.03394 

0.09006 

0.15884 

0.01860 

0.07070 

0.34948 

0.44391 

0.66305 

0.69091 

0.80851 
0 .89886 

0.  OfO.'S 


Circular  NACF.t.U 


Ci  ivumtVrrn  tial  Panel 
Rises  in  Terms  of  the 
An^le  0 in  Decrees 

0.0 

18.0 

36 .0 

64.0 

72.0 

90.0 

108.0 

1 06 .0 

144.0 

160.0 

180.0 


TABLE  3 

Presentation  of  datumresults 


B)  Data  key 

The  datum  results  are  given  in  part  C of  this  table.  The  results  pertaining  to  some  test  case  are  sub- 
divided in  one  or  more  of  the  following  data-blocks. 

- Case-identification  bloc 
RUBBERT 
DATUM 
RAE  WING 

.02 T/C  (or  C/d_  for  NACELLE  cases) 

5.0 ALFA 


- Pressure— data  block 

CP-DATA 

3_ 



_ n$  number  of 

40  - 

— 

- 

. mj  " " 

.079 

.549 

.924 

”i 

.9937  .04855 

.05068 

.05394 

X/C  i CPCn.) 

.9690  .OOO65 

.00291 

.01226 

.9207  -.02209 

-.01908 

-.00421 

- Veloci tv-data  block 

• 

• 

chordvd.se  ” ,f 

i = l(l)n 
i = l(l)n 


VELOCITY-DATA 

1 _ _ _ _ nj  number  of  spanwise  stations  for  which  data  are  given 

40 mj  " " chordwise  " " " " " " 

•549  \ 

.9937  .92978  .00269  -.09288 

.9690  .98301  -.01936  -.05578  x/c  , V^n.),  v/r,.),  V2(ti.)  i = i(i)n 

.9207  . 99988  -.02663  -.05373 

! j • i 

<Dragloog^data_jDl<5ck  i , 

DRAGLOOP-DATA 

1 — n,  number  of  spanwise  stations  for  which  data  Eire  given 

40 m,  " " chordwise  " " " " " " 

.549  ^ i = l(l)n 

.0006  .12688 

.0024  .03020  z/c  f cp(r,.)  i = l(l)n 

.0049  -.00336 

t • 

• I 


1 

ROBERT1.  , 

RAE  WING 

9 

t/c 

.15 

9 

ALFA 

.0  (CP-,  VELOCITY-DATA) 

2 

N 

9 

H If 

t 

t/c 

- 

.15 

9 

ALFA 

5.0  (CP-,  DRAGLOOP-DATA) 

3 

it 

9 

It  N 

9 

t/c 

- 

.05 

9 

ALFA 

= 

5.0  (CP-,  VELOCITY-,  DRAGLOOP-,  GAMMA-DATA) 

4 

•1 

9 

I*  II 

9 

t/c 

- 

.02 

9 

ALFA 

= 

5.0  (CP-,  VELOCITY-,  DRAGLOOP-,  GAMMA-DATA) 

5 

it 

9 

ST  RAKED  WING 

9 

t/c 

.05 

9 

ALFA 

5.0  (CP-,  VELOCITY-,  GAMMA-DATA) 

6 

It 

9 

it  t« 

9 

t/c 

= 

.02 

9 

ALFA 

= 

5.0  (CP-,  VELOCITY-,  GAMMA-DATA) 

7 

n 

9 

NACELiaE 

9 

c/de 

“ 

1.0 

9 

ALFA 

= 

.0  (CP-DATA) 

8 

it 

9 

it 

9 

c/de 

= 

3.333 

9 

ALFA 

= 

.0  (CP-DATA) 

9 

RUBBERT, 

RAE  WING 

9 

t/c 

= 

.05 

9 

ALFA 

= 

5.0  (CP-,  VELOCITY-,  DRAGLOOP-,  GAMMA-,  LOAI>-DATA) 

10 

N 

9 

11  it 

9 

t/c 

- 

.02 

9 

ALFA 

= 

5.0  (CP-,  VELOCITY-,  DRAGLOOP-,  GAMMA-,  LOAD-DATA) 

11 

M 

9 

STRAKED  WING 

9 

t/c 

= 

.02 

9 

ALFA 

= 

5.0  (VELOCITY-,  CP-,  GAMMA-,  LOAD-DATA) 

12 

H 

9 

NACELLE 

9 

C/De 

= 

1.0 

9 

ALFA 

= 

.0  (CP-DATA) 

13 

It 

9 

1* 

9 

o/de 

= 

1.0 

9 

ALFA 

= 

5.0  (CP-DATA) 

14 

HESS  , 

11 

9 

0/de 

= 

1.0 

9 

ALFA 

= 

.0  (CP-DATA) 

15 

w 

9 

11 

9 

o/de 

■= 

3.333 

9 

ALFA 

= 

.0  (CP-DATA) 

TABLE  3 (continued) 
Presentation  of  datumreaults 
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Oamma-data  block 

GAMMA-DATA 

12 nj  number  of  spanwise  stations  for  which  data  are  given 

.0245  .14528 

.078  .14528  ti  , T 

.146  .14645 

Load^da£a_block  • 

LOAD-DATA 

12  - n»  number  of  spanwise  stations  for  which  data  are  given 

.0245  .14359 

.078  .14463  ri  j cllc/2 

.146  .1455  29 
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TABLE  3 

PRESENTATION  OF  DATUM  RESULTS  (CONTINUED) 


C)  DATUM  RESULTS 


«oar«Ts 

OATU*  1 

9AI  WIN6 
0.13 
0.0 

CP-DATA 


0.079 

0.599 

0.929 

1.00376 

0.35696 

0.36566 

0.37500 

0.99631 

0.25966 

0.26639 

0.27720 

0.9007? 

0.16263 

0.19651 

0.20025 

0.«»5-*26 

0.12900 

0.19039 

0.19295 

0. 90006 

0.02757 

0.05019 

0.06512 

O.A??9? 

.0.09216 

-0.02929 

•0.00721 

0.73093 

-0.12053 

-0.09170 

•0.07261 

0.63067 

-0.16170 

-0.15666 

•0.13019 

0.52772 

-0.25637 

-0.22309 

•0.20660 

0.92665 

-0.29026 

-0.26962 

-0.27730 

0.33229 

-0.33369 

-0.35259 

-0.39976 

0. 29607 

-0.36939 

-0.90916 

-0.90000 

0.17655 

-0.37592 

-0.93553 

-0.93990 

0.11665 

-0.36029 

-0.9  3fc95 

-0.93097 

0.07966 

-0.30966 

-0.39367 

-0.39909 

0.09353 

-0.20963 

-0.29936 

-0.30022 

0.02269 

-0.09095 

-0.11519 

-0.12252 

0.01055 

0.16990 

0.13059 

0.12517 

0.00909 

0.91657 

0.37693 

0.37353 

0.00112 

0.56953 

0.59516 

0.53657 

0.00010 

0.63311 

0.61236 

0.60562 

1.00576 

0.35696 

0.56566 

0.37560 

0.99631 

0.25966 

0.26639 

0.27720 

0.96072 

0.16263 

0.19651 

0.20025 

0.95926 

0.12900 

0.19039 

0.15295 

0.90006 

0.02757 

0.03019 

0.06512 

0.62292 

-0.05216 

-0.02929 

-0.00721 

0.73093 

-0.12053 

-0.09170 

-0.07201 

0.63067 

•0.16170 

-0.15666 

-0.13019 

0.52772 

-0.23837 

-0.22309 

-0,20660 

0.92665 

-0.29026 

-0.26962 

-0.27730 

0.33229 

-0.33369 

-0.35259 

-0,39976 

0.29607 

-0.36959 

-0.90916 

-0,90066 

0.17655 

-0.37592 

-0.93353 

-0,93990 

0.11685 

-0.36029 

-0.93693 

-0.93097 

0.07988 

-0.30986 

-0.39367 

-0.39909 

0.09353 

-0.20963 

-0.29936 

-0.30022 

0.02289 

-0.09095 

-0.11519 

-0.12252 

0.01055 

0.18990 

0.13039 

0.12517 

0.00909 

0.91657 

0.37693 

0.37353 

0.00112 

0.56953 

0.59316 

0.53637 

0.00010 

0.63311 

0.61236 

0.60562 

vruoc I TV-OAT  A 

1 

92 

0.599 

1.00576 

0.76199 

0.07097 

-0,13291 

0.99631 

0.69277 

0.03332 

-0.19330 

0.96072 

0.66316 

0.03663 

-0.19073 

0.95926 

0.91999 

0.0239'. 

-0.15102 

0.90006 

0.96300 

0.00263 

-0.19972 

0.62292 

I .00176 

-0.01900 

-0.29320 

0.73093 

1.03567 

-0.03122 

-0.13503 

0.63067 

1.06600 

-0.09757 

-0.11029 

0.52772 

1 .09981 

•0.06963 

-0.09659 

0.92665 

1.13090 

•0.00202 

-0.06995 

0.33229 

1.15671 

-0.09013 

-0.01096 

0.29607 

1.17906 

-0.10900 

0.09306 

0.17655 

1.16691 

•0.11202 

0.12539 

0.11885 

1.17389 

-0.10195 

0.21916 

0.07966 

1.13266 

-0.07061 

0.32626 

0.09353 

1.09900 

-0.00001 

0.99056 

0.02269 

0.90999 

0.09993 

0.52751 

0.01055 

0.72336 

0.23269 

0.59057 

0.00909 

0.59003 

0.36702 

0.93950 

0.00112 

0.92171 

0.95319 

0.26600 

0.00010 

0.37390 

0.99073 

0.00593 

1.00576 

0.78199 

0.07097 

-0.13291 

0.99631 

0.69277 

0.03332 

-0.19330 

0.96072 

0.68318 

0.03663 

-0.19073 

0.95926 

0.91999 

0.02399 

-0.15102 

0.90006 

0.96300 

0.00263 

-0.19972 

0.6229? 

1.00176 

-0.01900 

-0.19326 

0.73093 

1.03537 

•0.03122 

-0.13505 

0.63067 

1.06600 

-0.09757 

-0.11629 

0.52772 

1.09981 

-0.06963 

-0.09659 

0.92665 

1.13090 

•0.00202 

-0.06995 

0.33229 

1.15671 

*0.09013 

-0.01096 

0.29607 

1.17906 

•0.10900 

0.09566 

0.17655 

1.18691 

•0.11202 

0.12359 

0.11665 

1.17369 

-0.10195 

0.21916 

0.07966 

1.13266 

•0.07061 

0.52026 

0.79353 

1.09900 

•0.00001 

0.99056 

0.02269 

0.90999 

0.09995 

0.52751 

0.01035 

0.72336 

0.23269 

0.59057 

0.00909 

0.59003 

0.36702 

0.93950 

0.00112 

0.92171 

0.95519 

0.26600 

0.06010 

0.37590 

0.99073 

0.06593 

.07966 

-1.00101 

.09355 

•1.15901 

.02269 

-1.13993 

.01055 

-.97772 

.00909 

-.65991 

.00112 

-.29900 

.00010 

.01903 

1.00576 

.37309 

.99651 

.2709? 

.96072 

.21996 

.95926 

.16669 

.90000 

.09211 

.02292 

.03505 

.73093 

-.01016 

.63067 

-.09075 

.52772 

-.06265 

.92665 

-.10935 

.33229 

-.12260 

.29607 

-.11907 

.17655 

-.07931 

.11005 

.00639 

.07900 

.12606 

,09353 

.29529 

.02269 

.97792 

.01055 

.60599 

.00909 

.39916 

.00112 

.93793 

.00010 

.25306 

OftAGLDOP-OATA 

9? 

0.599 

0.000 56 

0.36591 

0.00222 

0.26910 

0.00992 

0.16063 

0.00055 

0.12539 

0.01606 

0.02021 

0.02969 

-0.07256 

0.09196 

-0.16210 

0.05377 

-0.23370 

0.06369 

•0.35796 

0.07117 

-0.96976 

0.07979 

-0.59121 

0.07920 

-0.71005 

0.06976 

-0.09696 

0.06215 

•0.97216 

0.05297 

-1.00101 

0.09166 

•1.15401 

0.03139 
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TABLE  4 

Computers  used  and  their  relative  speeds  (estimated) 


COMPUTER 

USED 

SCALE  FACTOR  USED 

TO  OBTAIN  CDC  6600  TIME 

RUBBEHT 

CDC  6600 

1.0 

ROBERTS 

ICL  1906S 

.6 

HUNT  - SEMPLE 

IBM  370/158 

1.5 

NLR 

CDC  Cyber  7214 

3.0 

(e.g.  NLR  : 1 hr  CDC  6600  - 3.0  hrs  Cyber  7214) 


TABLE  5A 

Comparison  of  calculation  times 


METHOD 

PANEL  DISTRIBUTICH 

ON  HALF  NINO 

(W-'TOITC 

(FIGURES  HFHVEFH  BRACKETS  ARE  CDC  - 

time;  according  to  table  4) 

6600 

CHORD 

SPAN 

HAKE 

TIP 

NUMB.  OF 

MODES 

TOT.  NUMBER  OF 

SINGULARITIES 

AF, ROBYN  . 

INFL.  COEFF . 

SOL.  OF  UN 

SYST.  OF  ESS 

NUMBER  OF 

ITERATIONS 

TOTAL 

ROBERTA  1 2 WTUM 

3*1 

13 

13 

630 

11.90(39.74) 

8.9  (14.80) 

39.6  (117.33) 

27 

9 

9 

330 

4.95(16.5  ) 

1.35(  2.30) 

15.6  ( 49.33) 

19 

6 

6 

176 

2.60(  8.70) 

.3  ( .5  ) 

8.5  ( 28 .00) 

HUBRKRT(DATOM) 

40 

12 

12 

40 

684 

6.36(  6.36) 

2.53(  2.53) 

U.36(  11.36) 

22 

8 

8 

22 

296 

1 .37(  1.37) 

.32(  .32) 

2.46(  2.46) 

12 

12 

12 

12 

264 

1.00(  1.00) 

•24(  .24) 

1.90(  1.90) 

HUNT  - SKMTI.E 

60 

12 

12 

732 

31.2  (20.80) 

8.0  ( 5.33^ 

54 

40.0  ( 26.67) 

(SHEETS) 

30 

12 

12 

372 

8.47(  5.65) 

2.17(  1.45T 

56 

10.9  ( 7.27) 

NLR 

90 

12 

12 

1092 

48.94(16.31) 

4.6o(  1.53) 

12 

60.0  ( 20.0) 

60 

12 

12 

732 

21.66(  7.22) 

2.44(  .81) 

14 

28.38(  9.46) 

30 

12 

12 

372 

6.02(  2.00) 

•72(  .24) 

15 

Q.06(  3.0?) 

1)  Roberts  took  advantage  of  the  non-  cambered  nature  of  the  wings,  and  the  axi-symmetry  of  the  nacelle. 
This  is  accounted  for  in  the  figures  presented  between  brackets.  The  other  participants  employed  only 
the  XZ-plane  as  plane  of  symmetry, 

2)  Por  the  RAE  WING  with  T /C  » .05  these  times  are  less  than  one  half  those  for  T/C  - .02 


CASK:  RAE  KINO 
T/C  . .02 
a -5.0° 


-2>- 


|CASE:  STRiKED  WING 
T/C  - .02 
a - 5.0 


TABLE  5B 

Comparison  of  calculation  times. 


1 

CPU  - MINUTES 

PANEL  DISTRIBUTION 

(FIGURES  BETWEEN  BRACKETTS 

ARE  CDC  - 

6600  TIMES 

METHOD 

CN  HAU  WING 

ACCORDING  TO  TABLE  4) 

NUMB  .OF 

TOTAL  NUMB. OF 

AERODYN 

SOL. OF  LIN. 

NUMB.  OF 

CHORD 

SPAN 

WAKE 

TIP 

MODES 

SINGULARITIES 

INFL.COEFF . 

SYST.OF  BQS. 

ITERATIONS 

TOTAL 

ROBERTS 

39 

18 

966 

24.8  (82.67) 

31.02(51.73) 

81.3  (219.33) 

HUBBERT 

38 

12 

12 

38 

734 

6.91(  6.91) 

3.17(  3.17) 

12.55(  12.55) 

(DATOM) 

22 

8 

8 

22 

371 

1.96(  1.96] 

.55 ( -55) 

3.43(  3.43) 

12 

12 

12 

12 

298 

1.20(  1.20] 

.33(  .33) 

2.19(  2.19) 

HUNT  - SEMPLE 

60 

12 

12 

732 

31.30(20.87) 

8.8o(  5 .87 9 

6(J> 

40.8  ( 27.20) 

(SHEETS) 

30 

12 

12 

372 

8.83(  5.89) 

2.32(  1.55 9 

11.39(  7.59) 

NLR 

90 

12 

12 

1092 

49.20(16.40] 

4.15(  1.38) 

11 

59-83(  19.94) 

60 

18 

18 

1046 

42.94(14.31] 

7.51(  2.54) 

2C^ 

55-69(  18.56) 

60 

12 

12 

732 

22.17(  7.39] 

2.44(  .81) 

14 

28.93(  9.64) 

30 

12 

12 

372 

6.13(  2.04] 

.67(  .22) 

14 

9.17(  3.06) 

1)  90  solution  not  fully  converged. 

2)  Not  fully  converged. 

3)  For  T/C  = .05  these  times  are  less  than  one  half  those  for  T/c  = .02. 


table  5c 

Comparison  of  calculation  times 


CASE:  NACELLE 


1.0 


METHOD 

PANEL  DISTRIBUTION 

ON  HALF  NACELLE 

CPU  - MINUTES 

(FIGURES  BETWEEN  BRACKETS  ARE  CDC  - 6600  TIMES 
ACCORDING  TO  TABLE  4) 

CHORD 

CIRCUMF. 

WAKE 

NUMB  .OF 

MODES 

TOTAL  NUMB  .OF 

SINGULARITIES 

AERODYN 

INFL.COEFF. 

SOL. OF  SYST. 

OF  LIN.EQS 

NUMBER 

OF  ITER. 

TOTAL 

TIME 

ROBERTS( DATUM) 

55 

10 

58 

1.10(14.67) 

.033(  .053) 

2.90(38.0  ) 

RUBBERT(DATUM) 

40 

10 

10 

516 

3.74(  3.74) 

1.17  (1.17  ) 

6.52(  6.52) 

HUNT  - SEMPLE 

(sheets) 

60 

10 

10 

610 

21.50(14.33) 

2.60  (1.93  ) 

25 

24.6  (16.4  ) 

NLR 

60 

10 

10 

610  ||l4.96(  4.99) 

1.59  ( .53  ) 

15 

19.35(  6.63) 

• AM  BOUNDA»> 
CONDITION  POINTS 
SUM  AM  BOUNDAM' 
CONDITION  POINTS 


SUPTACf  SOURCE  PANELS 


F.g.  1 Def.m  t»on  of  (low  region  (or  3 D lifting  (low  probli 


HORSE  SHOE 
VORTEX  STSTI 


Schematic  view  of  selected  system  of  singularity 
distributions  for  the  NLR  panel  method 


PARABOLIC  SHAPE 


X PANEL  EDGES 


2 Chordwise  parabolic  shape  of  internal  vortex  strength 
variation  used  in  the  NLR  panel  method 


Sc  • normal  velocity  influence  coefficient  matrix  of 
sources  on  surface  collocation  points 

Sk  - refers  to  the  influence  of  sources  on  Kutta 
points 

Vc  t refers  to  the  influence  of  vortices  on  surface 
collocation  points 

Vk  refers  to  the  influence  of  vortices  on  Kutta 
points 

(T  unknown  source  singularities 

I =•  unknown  vortex  or  doublet  strength 

Ro  - boundary  condition  at  surface  collocation 
points 

Rk  boundary  condition  at  Kutta  points 


KUTTAPOINT 


X SURFACE  POINTS 
Q CAMBER  POINTS 


Definition  of  kutta* point  location  in  the  NLR  panel 
method  (N.B.  the  trailing  edge  is  open) 


Structure  of  normal  velocity  influence  coefficient  matrix, 
unknown  singularity  vector  and  right-hand  side  vector  in 
the  NLR  panel  method 
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FIG.  73  COMPARISON  OF  DATUM  RESULTS 
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FIG.  16  COMPARISON  OF  DATUM  RESULTS 
Y-COMPONENT  OF  VELOCITY 
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rlG.  2V  COMPARISON  OF  DATUM  RESULTS 
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FIG.  Z3  COMPARISON  OF  DATUM  RESULTS 

LHORDWISE  PRESSURE  DISTRIBUTION 

STRAKED  WING  . T/C  = .02  . a = 5.0  . n - .219 


( 


COMPARISON  OF  DATUM  RESULTS 
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FIG.  5/  FIRST  ORDER  METHOD  COMPRR1SON 
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FIG.  97  COMPARISON  OF  DATUM  RESULTS 
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